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Recently,  a  numerical  simulation  code  was  developed  to  study  the  interaction  of  the  effects  of  a  row  of 
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the  contribution  of  the  small  scales  by  a  subgrid  eddy  viscosity  model.  Using  both  coherent  (axisymmetric 
and  azimuthal)  and  random  excitation  of  the  impinging  jet,  the  evolution  of  the  large-scale  structures  in  the 
flow  Held  was  investigated  and  the  results  were  shown  to  agree  quite  well  with  experimental  data. 

In  il.e  present  Phase  I  study,  this  numerical  code  was  extended  to  include  other  realistic  V/STOL 
configurations,  such  as  an  obliquely  impinging  jet  and  an  impinging  jet  in  crossflow  with  and  without  a 
moving  ground  plane.  These  configurations  mode!  the  V/STOL  aircraft  in  pitch  and  in  forward  motion 
These  modifications  were  successfully  achieved  and  preliminary  simulations  using  axisymmetric  forcing  of  the 
impinging  jet  were  carried  out  to  demonstrate  the  capability  of  the  numerical  code.  Even  with  the  coarse 
mesh  (and  thus  low  Reynolds  number)  simulations,  the  results  show  good  qualitative  agreement  with 
experimental  observations.  For  example,  the  oblique  impinging  jet  simulation  showed  that  a  significant 
portion  of  the  jet  flow  is  carried  in  the  forward  direction  and  the  stagnation  point  on  the  ground  plane  is 
offset  from  the  jet  centerline.  The  shift  in  the  stagnation  point  from  its  location  for  a  normal  impinging  jet 
was  determined  to  agree  quite  well  with  experimental  observations,  and  the  presence  of  vorticity  in  the 
oblique  stagnation  streamline  was  demonstrated  which  also  agrees  with  past  observations.  The  upwash 
fountain  resulting  from  the  collision  of  the  wall  jets  on  the  ground  is  skewed  in  the  forward  direction  due  to 
the  obliqueness  of  the  impinging  jets.  This  would  result  in  a  nonuniform  pressure  on  the  undersurface  of  the 
V/STOL  aircraft  and  thus  cause  a  nonuniform  lift-off  effect. 

The  simulation  of  the  impinging  jet  in  crossflow  showed  the  formation  of  the  counterrotating  streamwisc  vor¬ 
tices,  which  has  been  observed  in  many  experiments.  However,  it  was  shown  that  the  presence  of  the  adja¬ 
cent  impinging  jet  (which  is  absent  i q  all  experiments)  can  significantly  modify  the  dynamics  of  the  vortical 
motion  in  the  flow  field.  Significant /three-dimensional  vortex  stretching  occurs  in  this  flow  and  the  redistri¬ 
bution  of  the  total  vorticity  in  the  flow  field  occurs  in  a  complex  manner.  The  simulation  with  the  ground 
plane  moving  at  the  same  velocity  as  the  crossflow  results  in  the  elimination  of  the  wall  boundary  layer.  As  a 
consequence,  the  size  and  the  structure  of  the  vortex  ahead  of  the  impinging  jet  is  drastically  reduced,  which 
is  similar  to  experimental  observations. 

The  results  presented  in  this  Phase  I  study  indicate  that  the  primary  technical  objectives  of  the  study  were 
successfully  achieved  and  a  numerical  simulation  capability  now  exists  that  can  be  utilized  to  understand  the 
dynamics  of  a  variety  of  flight  configurations  of  the  V/STOL  aircraft  in  ground  effect  in  a  systematic  manner. 
This  simulation  code  can  complement  experimental  studies  of  such  flow  fields  and  may  help  in  optimizing  the 
design  and  performance  of  V/STOL  aircraft. 
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1.  INTRODUCTION 


The  fluid  dynamics  of  impinging  jets  associated  with  vertical/short  takeoff  and  landing  (V/STOL) 
aircraft  operating  near  the  ground  is  very  complex.  An  impinging  jet  flow  field  is  a  combination  of 
free  jet  flow,  vortical  flow  in  curved  shear  layers,  stagnation  point  flow,  and  wall  jet  flow.  In 
addition,  the  effects  of  compressibility,  hot  gas  exhaust,  and  heat  transfer,  as  well  as  complex 
interactions  between  the  airframe  and  the  ground,  can  complicate  the  flow  field.  The  complexity  of 
the  flow  field  is  further  compounded  when  the  aircraft  is  in  forward  motion  near  the  ground  and  also 
when  the  aircraft  executes  a  pitch-up  or  pitch-down  motion  during  transition  in  and  out  of  ground 
effect.  Clearly,  an  accurate  prediction  of  the  aerodynamic  forces  and  moments  on  the  aircraft  is  not 
possible  without  understanding  the  associated  flow  physics. 

The  flow  fields  associated  with  the  various  operating  modes  of  V/STOL  aircraft,  such  as  hovering 
and  transition  in  and  out  of  ground  effect,  are  quite  different  from  each  other.  Many  of  the  complex 
phenomena  associated  with  these  How  fields  are  poorly  understood,  thereby  restricting  our  ability  to 
optimize  the  aircraft  design.  For  example,  lifting  jets  entrain  air,  w’hich  leads  to  induced  suction 
pressures  on  the  aircraft  undersurface.  When  the  aircraft  is  hovering  near  the  ground,  further 
entrainment  is  caused  by  the  wall  jets.  This  significantly  increases  the  suckdown  force  on  the  aircraft. 
The  problem  becomes  more  complicated  when  there  are  multiple  jets  interacting  near  the  ground,  as 
with  the  exhausts  of  a  twin-engine  V/STOL  aircraft.  The  wall  jets  caused  by  the  impinging  jets 
collide  and  form  an  upwash  or  a  fountain  that  impinges  on  the  aircraft  undersurface.  This  upwash 
creates  an  upload  that  partially  offsets  the  suckdown  force.  In  actual  flight  conditions,  however,  the 
upwash  of  the  hot  jet  exhaust  will  create  other  problems  such  as  hot  gas  ingestion  by  the  engine 
intakes  and  aircraft  skin  damage  caused  by  the  corrosive  effect  of  the  hot  combustion  products  in  the 
exhaust  and  by  the  entrainment  of  ground  surface  particles  by  the  wall  jets  In  general,  the  complex 
flow  field  associated  with  multiple  jets  in  ground  effect  is  not  very  well  understood. 

Most  of  the  research  that  has  been  directed  toward  understanding  V/STOL-related  flow  fields  has 
been  experimental. 1-5  Such  studies,  however,  addressed  mostly  global  features  and  time-averaged 
measurements.  Measurements  in  the  flow  field  associated  with  impinging  jets  is  extremely  difficult 
due  to  the  high  levels  of  turbulence  in  the  globally  unsteady  three-dimensional  flows.  Characterization 
of  these  flow's  would  require  measurements  using  three-dimensional  arrays  of  sensors,  such  as  hot  w  ire 
sensors,  which  are  limited  by  the  fact  that  they  cannot  discriminate  flow  reversal  from  forward  flow 
and/or  are  prone  to  probe  interference.  Thus,  it  is  clear  that  measurements  in  the  flow  field  of  three- 
dimensional  impinging  jets  will  require  new  measurement  techniques  such  as  double-pulse  holograph} 
or  scanning  two-color  LDA,  which  are  currently  being  developed. 

Numerical  simulation  provides  an  opportunity  to  study  the  detailed  flow  physics  as  a  function  of 
space  and  time.  Although  the  complete  flow  field  around  a  V/STOL  aircraft  cannot  yet  be  resolved 
adequately,  state-of-the-art  supercomputers  can  be  used  to  study  some  local  flow  domains  such  as  the 
impinging  jet.  In  addition,  numerical  simulation  has  the  additional  advantage  of  providing  the 
instantaneous  distribution  of  all  the  flow  variables  over  the  entire  computational  domain.  Thus,  it 
allows  "measurements"  of  flow  properties  not  possible  experimentally,  and  can  provide  simultaneous 
"flow  visualization"  and  "measurements"  in  arbitrary  planes.  Furthermore,  numerical  simulation  allows 
independent  control  of  flow  parameters  such  as  initial  condition,  free  stream  turbulence,  and 
excitation  frequency  and  amplitude.  These  parameters  can  be  more  easily  introduced  and  controlled 
in  a  numerical  simulation  than  in  laboratory  experiments. 
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In  this  Phase  I  study,  we  investigated  the  configurations  of  a  V/STOL  aircraft  in  pitch  while 
hovering  and  in  forward  motion  during  takeoff  and  transition.  The  dynamics  of  the  flow  in  the  foun¬ 
tain  caused  by  the  collision  of  the  wall  jets  will  be  affected  by  the  orientation  of  the  impinging  jet  and 
also  by  the  forward  motion  of  the  aircraft.  The  consequent  changes  in  the  lift-off  effect  on  the 
V,  STOL  aircraft  can  be  quite  significant.  A  detailed  understanding  of  the  dynamics  of  the  flow  field 
is  essential  to  optimize  the  design  of  the  aircraft.  To  study  the  complex  flow  field  associated  with  this 
configuration  in  ground  effect,  a  numerical  simulation  technique  must  be  developed  that  can  be  used 
to  understand  the  physics  of  the  flow.  The  demonstration  of  such  a  numerical  predictive  capability 
was  the  primary  objective  of  this  Phase  I  study. 

In  the  past,  extensive  experimental  and  numerical  studies  of  normal  impinging  jets  have  been  car¬ 
ried  out.3-5-10  Studies  of  inclined  impinging  jets  have  also  been  carried  out  in  the  past.5,  11  —  13  The 
earliest  work  by  Taylor1415  was  an  analytical  and  experimental  study  of  an  oblique  impinging  jet.  He 
showed  that  for  two-dimensional  impinging  jets,  the  maximum  (stagnation)  pressure  at  the  impinge¬ 
ment  point  is  the  same  for  all  angles  of  incidence  but  that  the  area  over  which  the  high  pressure  acts, 
as  well  as  the  downward  force  exerted  on  the  plate,  are  reduced  significantly  as  the  angle  becomes 
small.  He  also  developed  analytical  expressions  for  the  amount  of  mass  that  is  transported  in  the 
downstream  direction.  Later  experimental  study  by  Donaldson  and  Snedekcr5  of  normal  and  oblique 
impinging  jets  provided  data  on  the  mean  flow  characteristics  such  as  velocity  and  pressure.  They 
also  determined  the  azimuthal  distribution  of  the  wall  jet  and  showed  that  the  mass  and  momentum 
flux  is  concentrated  near  the  centcrplanes  but  is  skewed  depending  upon  the  angle  of  incidence.  The 
stagnation  streamline  is  also  offset  from  the  jet  physical  centerline  due  to  oblique  impingement. 

In  numerical  studies,  Rubel1617  investigated  the  oblique  impingement  of  jets  using  an  inviscid. 
rotational  flow  model  formulated  in  terms  of  vector  and  scalar  potentials  and  the  vector  vorticity.  He 
showed  that  the  oblique  impingement  of  a  round  jet  creates  a  fully  three-dimensional  flow  field  for 
which  a  strcamfunction  formulation  is  inappropriate.  It  was  also  shown  that  whereas  normal  impinge¬ 
ment  of  a  jet  produces  a  stagnation  line  with  zero  vorticity,  the  oblique  impingement  of  the  jet  results 
in  a  shift  in  the  stagnation  line  and  as  a  consequence,  the  stagnation  streamline  becomes  rotational.  In 
an  earlier  study,  H-’ves18  showed  that  in  the  vicinity  of  a  nearly  axisymmetric  rotational  stagnation 
region,  the  vorticity  exhibits  an  inverse  square  root  variation  with  the  vertical  distance.  The  numeii- 
cal  computations  agreed  with  this  prediction. 

Numerical  and  experimental  studies  of  jets  in  crossflow  have  been  going  on  for  a  long  time1 4 19-25 
due  to  a  wide  variety  of  engineering  applications  for  this  type  of  flow.  Examples  of  practical  applica¬ 
tions  include  film  cooling  of  turbine  blades,  transverse  fuel  injection  in  combustors,  pollutant  dispersal 
from  smoke  stacks,  and  V/STOL  aircraft  in  ground  effect.  Most  of  the  above-mentioned  studies  con¬ 
cerned  the  effects  of  free  jets  in  crossflow,  although  their  relevance  to  V/STOL  applications  arc  quite 
clear.  Recently,  a  numerical  study  of  an  impinging  jet  in  crossflow  was  carried  out  using  a  QUICK 
scheme.26  However,  this  investigation  focused  on  a  very  weak  crossflow  and  thus  does  not  give  a 
complete  picture  of  the  flow  field.  This  is  due  to  the  fact  that  the  dynamics  of  the  jet  interaction  w  ith 
a  crossflow  depends  quite  significantly  on  the  ratio  of  the  jet  velocity  to  the  crossflow  velocity.4 

In  all  the  experimental  studies  mentioned  above,  the  crossflow  toward  the  jet  was  generated  in  a 
wind  tunnel  with  a  fixed  ground  plane.  However,  to  model  a  realistic  V/STOL  aircraft  in  forward 
motion,  the  ground  plane  should  also  be  in  relative  motion.  Thus,  the  effect  ~r  a  moving  ground  plane 
must  be  understood  before  the  experimental  data  can  be  correlated  with  actual  flight  data.  Recent 
experiments27'28  have  studied  the  effect  of  a  moving  ground  plane  either  by  using  a  moving  ground 
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belt  or  by  moving  the  impinging  jet  itself.  These  studies  clearly  showed  that  there  are  significant 
differences  between  the  flow  fields  around  an  impinging  jet  when  the  ground  is  fixed  and  when  it  is 
moving.  For  example,  it  was  shown  that  the  vortex  strength  ahead  of  the  jet  is  reduced  by  as  much  as 
50  percent  in  lateral  penetration  when  the  ground  is  moving,  due  to  the  elimination  of  the  wall  boun¬ 
dary  layer.  In  addition,  actual  flight  test  data  indicates  that  there  is  less  lift  increase  during  approach 
and  landing  than  that  determined  from  wind-tunnel  measurements  with  a  fixed  ground  plane.  To 
evaluate  these  differences  experimentally  is  quite  difficult  and  expensive  since  special  test  facilities 
arc  required  to  simulate  the  moving  ground  plane.  However,  studying  the  effect  of  motion  of  the 
ground  plane  computationally  is  a  relatively  simple  matter  and  can  be  accomplished  through  appropri¬ 
ate  modifications  to  the  wall  boundary  conditions.  Numerical  investigation  of  this  configuration  of  a 
V/STOL  aircraft  in  forw-ard  motion  has  not  been  carried  out  before.  The  results  presented  in  this 
report  appear  to  be  the  first  attempt  to  model  the  effect  of  a  moving  ground  plane. 

In  this  Phase  I  study,  a  numerical  investigation  of  oblique  impinging  jets  and  impinging  jets  in 
crossflow  with  and  without  a  moving  ground  plane  has  been  carried  out.  This  involved  the  modifica¬ 
tion  of  a  three-dimensional  Navicr-Stokes  solver  that  was  recently  developed  to  study  impinging 
jets.0  10  The  problem  considered  in  our  previous  investigations  contained  the  essential  features  of 
twin  jets  impinging  normally  on  the  ground,  simulating  the  hovering  configuration  of  a  V/STOL  air¬ 
craft.  The  computer  code  solves  the  unsteady  three-dimensional  incompressible  Navicr-Stokes  equa¬ 
tions  using  the  technique  of  large-eddy  simulation  (LES).  The  technique  of  LES  is  well  known  and 
has  been  used  to  study  various  types  of  flow  fields  in  the  past.  In  terms  of  computational  accuracy  , 
this  technique  lies  between  the  Reynolds-Averaged  Navier-Stokes  (RANS)  and  the  direct  numerical 
simulations  (DNS)  techniques.  The  RANS  technique  essentially  solves  the  time-averaged  Navicr- 
Stokes  equations  using  turbulence  modeling  to  include  the  effect  of  turbulent  fluctuations.  This 
approach  docs  not  take  into  account  the  wide  variation  in  the  length  scales  in  the  turbulent  flow  field 
and  only  a  global  steady-state  solution  can  be  obtained.  In  contrast,  the  DNS  technique  attempts  to 
resolve  all  the  length  scales  in  the  flow  field  and  uses  no  modeling.  Temporal  as  well  as  spatial  accu¬ 
racy  is  maintained  in  the  DNS  simulations.  Even  with  state-of-the-art  supercomputers,  however,  only 
a  few  orders  of  length  scales  can  be  resolved  and  thus  the  DNS  technique  is  limited  to  v  ery  low  Rey¬ 
nolds  numbers.  On  the  other  hand,  the  technique  of  LES  resolves  all  the  length  scales  above  the  grid 
resolution  using  a  time-  and  space-accurate  numerical  scheme,  and  the  contribution  of  scales  smaller 
than  the  grid  resolution  arc  modeled  by  prescribing  a  sub-grid  eddy  viscosity  model  that  characterizes 
the  energy  cascade  into  the  unresolvable  scales.  In  principle,  given  a  reasonable  grid  resolution,  rela¬ 
tively  high  Reynolds  number  flows  can  be  studied  using  the  LES  technique. 

The  application  of  the  LES  technique  to  V/STOL-related  flow  fields  as  described  in  the  previous 
investigations910  was  the  first  attempt  to  use  LES  to  study  this  type  of  complex  flow  field.  In  these 
studies,  the  dynamics  of  the  evolution  and  the  propagation  of  large  coherent  structures  in  the  imping¬ 
ing  jet  shear  layer  and  in  the  wall-bounded  jet  was  studied  in  detail.  The  effect  of  introducing  a  dis¬ 
turbance  at  the  jet  exit  plane  on  the  evolution  of  the  coherent  structures  was  determined.  The  effects 
of  axisymmetric,  azimuthal,  and  random  perturbation  were  investigated.9,10'29  Good  agreement  w  ith 
experimental  data  was  obtained,  as  described  earlier.9'10 

In  this  study,  some  modifications  to  the  basic  LES  code  were  necessary  to  model  an  obliquely  imp¬ 
inging  jet  and  an  impinging  jet  in  crossflow  with  and  without  a  moving  ground  plane.  These  modifi¬ 
cations  were  successfully  accomplished  and  some  numerical  simulations  were  carried  out  to  demon¬ 
strate  the  capability  of  the  modified  LES  code.  This  report  documents  the  findings  of  this  study. 
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2.  PHASE  I  TECHNICAL  OBJECTIVES 

The  primary  technical  objectives  of  this  Phase  I  study  are  given  below: 

.  Modify  the  existing  three-dimensional  code  to  numerically  simulate  inclined  impinging  jets. 

.  Modify  the  code  to  numerically  simulate  impinging  jets  in  a  crossflow  with  and  without  a  mo\ 
ing  ground  plane. 

•  Demonstrate  the  capability  of  the  numerical  code  to  simulate  inclined  impinging  jets  and  imp 
inging  jets  in  a  crossflow. 

As  discussed  in  the  following  sections,  these  Phase  I  objectives  were  successfully  achieved. 
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3.  NUMERICAL  MODEL 


The  numerical  code  used  in  this  study  was  originally  designed  to  solve  the  problem  of  a  row  of  jets 
impinging  vertically  from  a  horizontal  plate  onto  a  ground  plane.  Before  presenting  the  details  of  the 
modifications  required  for  this  Phase  1  study,  we  will  briefly  describe  the  basic  numerical  scheme. 

3.1  Governing  Equations 

The  governing  equations  of  the  numerical  scheme  are  the  three-dimensional,  incompressible 
Navier-Stokcs  equations  and  the  continuity  equation.  Following  classical  LES  approach,  the  How 
\ariablcs  are  decomposed  into  a  large-scale,  resolvable  part  and  a  subgrid-scale,  unresolvable  part. 
The  equation  governing  the  large-scale  motion  is  obtained  by  applying  a  cell-volume  averaging  to  the 
Navicr-Stokcs  equations  and  to  the  continuity  equation.  The  resulting  equations  are 

Qt  +  ^\QQT  ~  2(i/  +  v.)  f]  =  —  VP  (1) 

and 

V-7  =  0  (2) 

where  ~q  is  the  resolvable  part  of  the  velocity,  5  is  the  rate-of-strain  tensor,  and  v  is  the  kinematic 
v  iscosity.  Here,  the  subscript  t  denotes  differentiation  with  respect  to  time.  The  continuity  equation 
[Equation  (2)]  is  then  replaced  by  a  Poisson’s  equation  for  the  "modified"  pressure  P,  which  now  con¬ 
tains  the  terms  appearing  due  to  the  volume  averaging: 

V2/>  =  -  V-ft  -  VQ  (3) 

where 

Q  =  V[qqT-2(v+v,)S]  (4) 

and 

V-qt-  0  (5) 

The  system  of  Equations  (I)  and  (4)  is  equivalent  to  the  original  system  of  Equations  (1)  and  (2)  and 
is  used  here  in  place  of  the  original  system. 

It  is  assumed  that  the  small-scale  structures  are  universal  in  character,  contain  little  energy,  and  arc 
dissipative  in  nature.  A  subgrid-scale  eddy  viscosity,  is  defined  to  model  the  cascade  of  energy 
from  the  large  scale  to  the  small  scale,  where  it  is  dissipated.  The  form  of  the  eddy  viscosity  is  similar 
to  the  model  used  by  Smagorinsky30  and  Deardorff31  and  is  given  as 

i/,  =  (Ci A)2  (2 SaS^I2  (6) 

where  Cx  is  a  constant  and  A  is  the  characteristic  length  scale  of  the  smallest  resolvable  eddies.  In  the 
present  simulations,  Ci  =  0.1,  which  is  the  value  chosen  by  Deardorff.31 

The  physical  problem  investigated  in  these  simulations  is  an  infinite  row  of  impinging  jets.  Each 
jet  exits  from  a  circular  nozzle  of  diameter  D  in  a  horizontal  plate  that  is  parallel  to  the  ground  and  at 
a  distance  //above  it.  This  plate  models  the  undersurface  of  the  V/STOL  aircraft.  Figures  la  and  lb 
show  the  ♦op  and  :ide  views,  respectively,  of  the  normal  impinging  jet.  A  Cartesian  frame  of  refer¬ 
ence,  x  =  (x,r,:).  is  defined  as  shown,  with  the  velocities  in  each  coordinate  direction  defined  as 

q  -  (lt,\\w). 
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3.2  Boundary  Conditions 

The  filtered  Navier-Stokes  equations,  Equations  (1)  and  (4),  are  solved  subject  to  appropriate 
boundary  conditions.  By  virtue  of  the  assumption  of  an  infinite  row  of  impinging  jets,  periodic  boun¬ 
dary  conditions  are  imposed  at  the  boundaries  x  =  xt  and  x  =  xr  (see  Figure  1).  Even  when  a  high- 
resolution  grid  is  used,910  it  is  still  not  fine  enough  to  resolve  the  viscous  sublayer  adjacent  to  the  top 
and  bottom  walls.  Therefore,  the  solution  in  the  near  wall  region  is  patched  to  the  solution  at  the  first 
computational  cell  next  to  the  wall  by  using  Spalding’s  law  of  the  wall.32  The  vertical  velocity  com¬ 
ponent  (vv)  is  specified  as  zero  at  both  the  upper  and  lower  boundaries,  except  at  the  jet  exit  where  it 
is  specified  as  a  function  of  space  and  time.  The  other  two  velocity  components  (u,  v)  are  specified  as 
zero  at  the  jet  exit. 

If  there  is  a  crossflow  present,  the  boundary  conditions  at  the  outflow  boundaries  (i.e.,  y  =  yt  and 
v  =  y r)  must  be  modified.  For  oblique  impingement,  the  conditions  at  the  jet  exit  plane  also  must  be 
modified.  These  differences  are  highlighted  in  the  following  sections 

3.2.1  Normal  Impinging  Jet 

The  outflow  boundaries  at  y  =  y(  and  y  =  yr  are  artificial  boundaries  used  to  limit  the  computa¬ 
tional  domain  to  a  finite  region.  Since  the  boundary  conditions  at  these  locations  are  not  known, 
approximate  boundary  conditions  are  imposed.  Zero  derivatives  for  the  velocity  components  u  and  u 
are  imposed  at  these  boundaries  (i.e.,  du/dy  =0,  dw/dy  =  0)  and  the  derivative  dv/dy  was  deter¬ 
mined  from  the  continuity  equation.  To  solve  the  Poisson’s  equation  for  pressure,  Neumann  boundary, 
conditions  were  used  on  all  boundaries.  At  the  top  and  bottom  walls,  these  conditions  were  defined 
by  the  .v-momentum  equation.  At  the  outflow  boundaries,  dP/dy  =  0  was  imposed.  However,  in 
some  earlier  calculations,  this  weak  boundary  condition  allowed  a  continuous,  nonphysical  buildup  of 
the  inflow  velocity  in  certain  regions  of  the  boundary.  The  pressure  gradient  at  the  outflow  boun¬ 
daries  was  modified  to  limit  this  build  up.  Also,  the  existence  of  a  solution  to  the  Poisson’s  equation 
for  pressure  with  Neumann  boundary  conditions  requires  that  the  integral  over  the  computational 
domain  of  the  source  terms  in  that  equation  be  equal  to  the  integral  of  dP/dn  over  the  boundaries, 
where  n  is  the  outward  normal  at  the  boundaries.  To  satisfy  this  condition,  the  value  of  dP/dn  was 
uniformly  modified  at  the  outflow  boundaries. 

Note  that  these  modifications  to  the  pressure  at  the  outflow  are  not  required  when  there  is  a 
crossflow  present.  With  a  crossflow,  the  zero  derivative  for  the  pressure  mentioned  above  is  used 
without  any  condition  to  limit  the  inflow  velocity  since  there  is  no  reverse  flow  at  the  outflow  boun¬ 
daries. 

3.2.2  Modifications  to  Model  an  Inclined  Impinging  Jet 

To  simulate  an  obliquely  impinging  jet,  some  modifications  of  the  boundary  conditions  at  the  jet 
exit  plane  were  required.  In  our  earlier  study,  the  horizontal  velocity  components  u  and  v  were  both 
specified  as  zero  in  the  plane  of  the  jet  exit.  To  simulate  an  inclined  jet,  a  nonzero  velocity  com¬ 
ponent  v  in  the  y-direction  must  be  imposed.  The  required  modification  is  described  below. 
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First,  an  axial  velocity  distribution  at  the  circular  jet  exit  of  a  noninclined  jet  is  specified.  Wc 
assume  that  this  circular  jet  is  inclined  in  the  y-direction  such  that  the  upper  horizontal  computational 
boundary  will  intersect  the  jet  circular  cylinder  in  an  elliptical  cross  section  with  its  major  and  minor 
axes  in  the  y-  and  x-directions,  respectively.  Figure  2a  shows  a  schematic  of  this  configuration.  At 
the  computational  mesh  points,  which  fall  within  this  ellipse,  the  original  axial  velocity  component 
may  now  be  decomposed  into  its  y-  and  c-components.  This  velocity  distribution  is  then  specified  at 
the  jet  exit  plane  and  thus  simulates  the  case  of  an  inclined  impinging  jet. 

3.2.3  Modifications  to  Model  Normal  Impinging  Jets  in  Crossflow 

To  simulate  the  effect  of  crossflow,  further  modifications  to  the  code  were  required.  We  assumed 
the  crossflow  to  be  specified  in  the  positive  y-direction,  normal  to  the  plane  of  the  jet  axes.  To  simu¬ 
late  the  case  of  forward  motion  of  the  V/STOL  aircraft  in  the  negative  y-direction  at  a  given  vclocitv 
v  =  v 0,  a  positive  velocity  v0  in  the  y-direction  is  first  imposed  on  the  flow  and  the  boundary  so  that 
the  problem  can  be  solved  in  a  reference  frame  in  which  the  aircraft  is  fixed  (see  Figure  2b).  The 
governing  equation  for  this  case  remains  unchanged  except  for  the  case  in  which  the  aircralt  is 
accelerating  or  decelerating.  In  this  case,  there  is  an  additional  term,  dv0/dt ,  which  appears  in  the  y- 
momentum  equations.  While  the  case  of  an  accelerating  (or  decelerating)  aircraft  is  of  general  interest, 
we  did  not  go  into  the  details  of  these  unsteady  effects  in  this  Phase  I  study.  However,  the  code 
modification  was  carried  out  to  include  this  unsteady  term,  so  that  an  option  exists  to  study  forward 
acceleration  or  deceleration  while  in  forward  motion.  In  the  present  study,  this  option  was  used  to 
ramp  up  the  initial  transitional  period  in  which  the  crossflow  builds  from  an  initial  value  of  zero  to  its 
final  desired  value. 

3.2.4  Modifications  to  Model  Crossflow  with  Moving  Ground  Plane 

As  mentioned  earlier,  it  is  of  great  interest  from  an  experimental  point  of  view  to  compare  the 
results  of  impinging  jets  in  crossflow  with  a  moving  ground  plane  to  those  with  a  fixed  ground  plane. 
For  the  case  of  an  impinging  jet  with  crossflow  over  a  fixed  ground  plane,  the  modifications  described 
above  are  applicable.  For  the  case  with  a  moving  ground  plane,  in  the  reference  frame  of  a  fixed  air¬ 
craft,  the  ground  plane  will  also  move  in  the  positive  y-direction  with  the  velocity  v  =  v0  (sec 
Figure  2b).  This  modification  has  also  been  successfully  accomplished. 

Although  the  above  modifications  for  the  inclined  jet  and  the  crossflow  have  been  described 
separately,  the  boundary  conditions  for  these  configurations  can  be  prescribed  together  without  any 
problem.  Thus,  in  principle,  it  is  possible  to  model  a  V/STOL  aircraft  in  pitch  while  in  forward 
motion.  This  flight  situation  is  also  very  interesting  in  that  the  dynamics  of  the  flow  field  can  be 
modified  significantly  due  to  the  combined  effect  of  the  inclined  impinging  jet  and  the  forward 
motion.  These  aspects  of  the  flow  can  be  studied  in  more  detail  in  the  next  phase. 

3.3  Numerical  Method 

The  governing  equations  are  solved  using  a  staggered  computational  mesh  which  alleviates  many  of 
the  difficulties  associated  with  the  solution  of  incompressible  Navier-Stokes  equations  on  a  regular 
mesh.  A  brief  description  of  the  discretization  of  the  various  flow  variables  (i.e.,  P,  u,  v,  and  »»•)  on  a 
staggered  mesh  is  given  here. 
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The  main  computational  mesh  is  the  P  mesh.  The  pressure  is  defined  at  the  nodal  points 
(a,,  fj,  zk)  =  (/Ax,  j Ay,  kAz)  of  this  mesh,  which  is  uniform  in  the  three  coordinate  directions.  The 
u-,  v-,  and  w-velocity  components  are  defined  at  the  nodal  points  (x,  +  1/2,  zk),  (x„  y3  +  l/2,  :k)  and 
(a,,  \-j,  zk  +  1/2)  of  the  u,  v,  and  w  meshes,  respectively.  These  meshes  are  displaced  from  the  P  mesh 
by  half  of  one  mesh  spacing  in  the  x-,  y-,  and  z-directions,  respectively.  The  x-,  y-,  and  z-momentum 
equations  are  discretized  at  the  nodal  points  of  the  u,  v,  and  w  meshes,  respectively,  while  the  discre¬ 
tized  Poisson’s  equation  for  pressure  is  obtained  at  the  nodal  points  of  the  P  mesh  by  central  dif¬ 
ferencing  of  the  discretized  momentum  equations  at  the  neighboring  points. 

The  pressure  gradient  terms  and  the  diffusive  terms  of  the  momentum  equations  are  approximated 
by  central  differencing.  The  use  of  central  differencing  to  approximate  the  convective  terms  leads  to 
nonphysical  oscillations  at  relatively  low  Reynolds  numbers  and  to  divergence  of  the  solution  as  the 
Reynolds  number  is  increased.  To  avoid  this  problem,  the  convective  terms  of  the  momentum  equa¬ 
tions  are  discretized  using  a  three-dimensional  version  of  the  QUICK  scheme  first  developed  for 
one-dimensional  flows  by  Leonard33  and  then  extended  to  two  dimensions  by  Davis  and  Moore.3'1 
The  details  of  the  three-dimensional  QUICK  scheme  used  in  the  present  code  was  described  earlier0 
and  thus  will  not  be  repeated  here. 

To  advance  the  solution  from  time  tn  to  time  /n+1,  where  tn  =  nAt  and  At  is  the  incremental  time 
step,  the  Poisson’s  equation  for  pressure 

j?  _,n+l  _  (  _n 

vVn+1/2  =  —  +  i _ on+1/2 

At  At 

is  first  solved  to  evaluate  the  pressure  distribution  that  will  allow  the  continuity  equation  to  be  satis¬ 
fied  at  time  in+1.  Here, 


g-+w.i.g*_ig— 


and  V  is  the  central  difference  approximation  of  the  operator  V.  For  the  continuity  equation  to  be 
satisfied  at  time  /n+1,  it  is  necessary  to  set  the  first  term  on  the  right-hand  side  of  Equation  (7)  equal 
to  zero.  An  efficient  method  for  direct  solution  of  the  discrete  Poisson’s  equation  is  used.  Once  the 
pressure  field  is  determined,  an  Adams-Bashforth  time-stepping  scheme  is  used  to  calculate  the  velo¬ 
city  at  the  next  time  level.  Therefore, 

qn+1  =  qn -  At  [e"+1/2  +  ^/>B+1/2)  (9) 


The  finite  difference  scheme  used  in  these  simulations  is  second-order  accurate  in  both  space  and 
time. 
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4.  APPROACH  TO  NUMERICAL  SIMULATIONS 


In  this  study,  the  primary  technical  objectives  were  to  implement  the  modifications  to  simulate 
inclined  impinging  jets  and  jets  in  crossflow,  and  to  demonstrate  the  capability  of  the  numerical 
scheme  to  study  flow  fields  associated  with  such  V/STOL  flight  configurations.  Due  to  the  resource 
and  time  constraints  of  the  Phase  I  study,  only  representative  simulations  were  performed  to  demon¬ 
strate  the  code’s  capabilities.  As  shown  below,  however,  analysis  of  the  computed  solutions  indicates 
interesting  features  that  are  in  good  qualitative  agreement  with  past  experimental  observations. 

The  approach  used  in  this  study  was  to  investigate  the  motion  of  large-scale  coherent  structures  in 
the  flow  field.  There  are  various  reasons  for  using  this  approach.  Large-scale  coherent  structures 
have  been  extensively  observed  in  free  jets36,36  and  in  impinging  jets.3  These  structures  are  formed 
when  the  unstable  jet  shear  layer  rolls  up  into  vortices  that  pair  or  merge  during  their  downstream 
propagation  to  form  the  large  coherent  structures.  Due  to  three-dimensional  effects,  these  structures 
eventually  break  down  into  small-scale  turbulent  flow.  The  formation  of  these  structures  occurs  at  a 
preferred  frequency,  often  called  the  jet  preferred  mode,  which  lies  in  the  Strouhal  Number 
(St  =  fD/U)  range  of  0.2  to  O.8.37  It  has  become  well  known  in  recent  years  that  these  structures  play 
a  major  role  in  the  transport  of  momentum  and  in  noise  production.6  There  is  also  evidence  that 
these  structures  play  a  major  role  in  the  entrainment  process.  This  is  particularly  relevant  to  V/STOL 
aircraft  operating  near  the  ground  since  the  vortical  structures  may  entrain  dust  from  the  ground  that 
may  in  turn  be  ingested  into  the  engine  intakes,  possibly  causing  malfunctions.  Understanding  the 
dynamics  of  the  evolution  and  motion  of  large-scale  coherent  structures  in  impinging  jet  flows  is 
clearly  important  in  optimizing  the  design  of  V/STOL  aircraft. 

In  addition,  the  LES  technique  is  uniquely  suited  to  the  study  of  the  motion  of  coherent  structures 
since  these  structures  are  usually  much  larger  than  the  grid  resolution  and  thus  can  be  resolved  by  the 
numerical  scheme.  As  mentioned  earlier,  only  the  scales  smaller  than  the  grid  resolution  need  to  be 
modeled,  which  can  be  accomplished  by  using  a  sub-grid  model. 

To  study  the  motion  of  these  coherent  structures,  the  jet  inflow  velocity  profile  at  the  jet  exit  is 
perturbed  by  superposing  prescribed  disturbances.  The  resulting  motion  of  the  flow  is  then  calculated 
using  the  LES  code.  This  method  of  prescribing  external  excitation  was  very  successful  in  our  earlier 
studies  of  normal  impinging  jets.9,10,29  In  addition,  the  method  of  external  forcing  at  a  specified  fre¬ 
quency  suppresses  the  natural  instability  of  the  flow;  hence,  the  dynamics  of  the  coherent  structure 
are  controlled  by  the  forcing  frequency.  This  results  in  a  "clean"  numerical  experiment.  Externally 
excited  flows  have  been  studied  extensively  in  the  past,  both  in  free  jets35,38,39  and  in  impinging  jets.3 
to  understand  the  dynamics  of  large-scale  coherent  motion. 

4.1  Relevance  of  Grid  Resolution 

To  resolve  all  the  relevant  scales  of  motion  associated  with  high  Reynolds  number  flows,  a  very 
high  grid  resolution  is  required.  In  our  earlier  studies,9,10,29  for  example,  grid  resolutions  on  the  order 
of  64  x  128  x  32  were  used.  Computations  using  such  high  resolution  can  only  be  performed  on 
supercomputers  such  as  the  Cray.  The  numerical  code  has  been  developed  to  take  advantage  of  vec¬ 
tor  processing  architecture  and  thus  is  reasonably  fast  on  the  Cray.9,10  In  the  present  Phase  I  study, 
however,  the  primary  objectives  were  to  demonstrate  that  an  inclined  impinging  jet  and  an  impinging 
jet  in  crossflow  can  be  modeled  by  the  code.  Since  coarse  grid  computations  are  sufficient  to  demon¬ 
strate  this  capability,  simulations  on  a  supercomputer  were  not  carried  out  in  the  Phase  I  study. 
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Numerical  calculations  were  performed  on  the  in-house  MassComp  computer.  Due  to  the  mcmor\ 
and  speed  limitations  of  this  scalar  processing  computer,  only  coarse  grid  computations  were  per¬ 
formed.  Nevertheless,  it  was  desirable  that  the  results  of  these  simulations  be  useful  in  understanding 
the  flow  Held  associated  with  inclined  impinging  jets  and  impinging  jets  in  crossflow.  Also,  wherever 
possible,  the  potential  for  comparison  with  available  experimental  data  was  desirable  in  order  to  vali¬ 
date  the  calculations  and  to  demonstrate  the  strength  of  the  LES  technique.  Thus,  after  some  experi¬ 
mentation,  a  grid  of  24  x  64  x  24  was  chosen  as  an  optimum  compromise,  considering  the  computer 
speed  and  resource  limitations.  The  computational  domain  was  modeled  to  extend  from  -1.5  <  x  /D  < 
1.5,  and  -3  <  y  /D  <  3,  and  a  jet-to-ground  distance  (H /D)  of  1.5  was  simulated.  Since  the  configura¬ 
tions  studied  in  Phase  I  involved  oblique  impingement  in  the  y-z  plane  and  the  crossflow  was  also 
prescribed  in  the  y-z  plane,  the  computational  domain  in  the  x-z  plane  (i.e.,  the  plane  of  the  fountain) 
could  be  reduced  by  assuming  symmetry  conditions.  Thus,  at  x  =  0,  the  symmetry  boundary  condi¬ 
tions  are  used  and  the  computational  domain  in  the  x-c  plane  is  reduced  to  0  <  x  /D  <  1.5. 

In  our  earlier  studies,910  it  was  determined  that  the  evolution  of  large-scale  structures  can  be 
effected  by  the  jet-to-ground  distance.  For  example,  when  the  jet  exit  plane  is  close  to  the  ground, 
the  ground  plane  can  act  as  a  damping  mechanism  on  the  instability  waves  in  the  impinging  jet  shear 
layer.  The  computational  domain  in  the  horizontal  directions  (x,  y)  is  also  important.  For  example,  in 
the  x-direction,  the  domain  determines  the  spacing  between  adjacent  impinging  jets,  which  in  turn 
affects  the  dynamics  of  the  fountain.  The  effect  of  varying  these  parameters  can  be  considered  in  the 
next  phase  using  a  higher  grid  resolution. 

For  the  grid  resolution  used  in  this  Phase  I  study,  a  Reynolds  number  of  300  was  assumed  based 
on  the  jet  diameter  D  and  the  unperturbed  jet  exit  velocity  W ,  both  of  which  were  chosen  to  be  unity. 

4.2  Initialization 

Here,  we  describe  the  initialization  procedure  used  in  the  present  computations.  The  initialization 
for  the  steady-state  computations  and  the  coherently  forced  computations  are  described  separately. 

4.2.1  Steady-State  Profile 

The  computations  were  begun  with  no  flow  in  the  domain.  At  t  =  0,  a  jet  with  a  profile  of 

Wo(r)  =  -»T[l  —  (— H  (10) 

ri 

was  introduced  at  the  jet  exit.  Here,  W  is  the  unperturbed  axial  jet  velocity  magnitude,  which  is 
specified  as  1;  r}  =  D  / 2  is  the  radius  of  the  jet;  and  r  is  the  radial  coordinate  with  respect  to  the  jet 
center.  Past  studies  of  free  jets  indicate  that  the  initial  jet  velocity  profile  has  a  "top-hat"  shape.4'1 
However,  to  resolve  the  jet  shear  layer  for  the  "top-hat"  profile  requires  a  very  high  grid  resolution, 
which  was  not  possible,  as  mentioned  earlier.  The  profile  chosen  for  the  present  calculations  was  an 
optimum  compromise  for  the  grid  resolution  used. 
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4.2.2  Coherent  Forcing 

To  study  the  effect  of  external  excitation  on  the  dynamics  of  the  impinging  jets,  a  perturbation  w' 
is  added  to  the  steady-state  profile  given  in  Equation  (10).  The  perturbation  is  of  the  form 

w’=  —A  w0(r)sin(V>)  (11) 

where  ip  =  2-kI /T,  T  is  the  time  period  of  the  forcing  frequency,  and  A  is  the  amplitude  of  forcing 
Typically,  the  amplitude  was  varied  from  10  to  20  percent  of  the  mean  velocity.  This  is  similar  to 
experimental  studies  of  forced  impinging  jets3  and  to  our  earlier  study  of  normal  impinging  jets.  A 
forcing  frequency  was  chosen  to  correspond  to  a  Strouhal  number  (St  =  fD/U)  of  0.47.  Her,.. 
U  =  IF  -  1  is  the  reference  velocity  of  the  unperturbed  jet.  The  Strouhal  number  chosen  fjr  these 
simulations  is  in  the  jet  preferred  range.37  Furthermore,  the  effect  of  forcing  at  this  frequency  was 
studied  earlier  for  the  case  of  normal  impinging  jets.8'10  In  this  Phase  I  study,  symmetry  conditions 
were  imposed  at  x  =  0  in  order  to  adequately  resolve  the  large-scale  structures  with  the  available  grid 
resolution.  This  implies  that  only  axisymmetric  forcing  of  the  jet  could  be  studied.  To  study  the 
effect  of  azimuthal  or  random  excitation  of  the  impinging  jet,10  the  full  jet  has  to  be  modeled.  This 
requires  a  much  higher  grid  resolution  and  therefore  will  have  to  await  a  future  study. 
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5.  RESULTS  OF  SIMULATIONS 

In  this  section,  we  describe  in  detail  the  results  of  the  numerical  study  carried  out  in  this  investiga¬ 
tion.  This  section  is  divided  into  four  parts  dealing  with  the  four  flow  configurations  of  interest:  the 
oblique  impinging  jet,  the  normal  impinging  jet  without  crossflow,  the  normal  impinging  jet  in 
crossflow,  and  the  normal  impinging  jet  in  crossflow  with  a  moving  ground  plane. 

5.1  The  Oblique  Impinging  Jet 

The  computations  of  the  oblique  impinging  jet  were  begun  with  no  flow  in  the  computational 
domain  and  at  i  =  0.  The  steady-state  jet  velocity  as  prescribed  by  Equation  (10)  was  imposed  at  the 
jet  exit. 

A  starting  vortex  is  shed  by  the  jet  during  its  initial  propagation.  This  vortex  impinges  on  the 
ground  plane  and  propagates  along  the  ground.  For  example,  Figure  3  shows  a  three-dimensional  per¬ 
spective  of  the  total  vorticity  M  in  the  jet  impinging  at  an  angle  of  60°  (see  Figure  2).  Vorticity  higher 
than  the  level  shown  will  be  contained  inside  the  surface  shown.  Note  that  in  all  the  thrcc- 
Jiincnsionai  perspective  views  shown  in  this  report,  the  full  computational  domain  is  shown  to 
improve  the  clarity  of  presentation.  The  vorticity  in  the  full  jet  is  generated  by  forming  a  symmetric 
image  in  the  -1.5  <  x  /D  <  0  domain  from  the  computed  solution  in  the  0  <  x  /D  <  1.5  domain  by 
using  the  symmetry  conditions  at  .x  =  0.  Figure  3  shows  the  time  when  the  starting  vortex  ring  has 
impinged  on  the  ground  and  is  propagating  along  the  ground.  In  our  earlier  study  of  normal  impinging 
jets.0  it  was  determined  that  the  phase  speed  of  this  ring  on  the  ground  is  initially  the  same  in  all 
radial  directions,  but  as  the  vortex  ring  approaches  the  fountain  axes,  the  rings  from  the  two  adjacent 
wall  jets  collide  and  rise  into  the  fountain.  This  forms  the  upwash  that  eventually  impinges  on  the  top 
wall,  resulting  in  the  so-called  "lift-off"  effect.  The  motion  of  the  vortex  in  the  fountain  reduces  the 
phase  speed  of  the  vortex  ring  on  the  ground  plane  and,  subsequently,  the  vortex  ring  starts  to  tw  ist 
and  eventually  breaks  up  into  isolated  structures.  When  the  computation  was  continued  with  the 
steady-state  profile  as  given  by  Equation  (10),  the  flow  field  eventually  reached  a  steady  state. 

For  the  inclined  jet,  the  vortex  ring  does  not  propagate  uniformly  in  all  directions  on  the  ground 
plane.  Depending  upon  the  angle  of  impingement,  only  a  relatively  small  portion  of  the  fluid  will  be 
transported  in  the  backward  direction.  For  example,  Taylor15  showed  from  conservation  of  momen¬ 
tum  parallel  to  the  plate  that,  for  a  two-dimensional  jet  impinging  at  an  angle  9,  the  fraction  of  the 
fluid  that  flows  forward  is  cos20/2  and  the  fraction  that  flows  backward  is  sin20/2.  This  implies  that 
for  an  angle  of  60°,  only  25  percent  of  the  fluid  is  transported  in  the  backward  direction.  In  the 
three-dimensional  case,  this  is  not  true;  the  mass  flux  in  the  forward  and  the  backward  direction  can 
be  determined  only  by  integrating  the  total  flux  in  the  azimuthal  direction.  Rubel17  derived  an 

expression  for  the  fraction  of  fluid  carried  forward  in  the  three-dimensional  case  and  showed  that  this 
0  3 1  n  ^  0 

fraction  is  ( 1  — - 1-  — — ),  which  for  the  present  case  of  9  =  60°,  gives  a  value  of  around  0.805.  This 

7T  2n 

value  agrees  quite  well  with  the  experiments  ofSchach13  and  Taylor.14 

The  dynamics  of  the  vortical  structures  being  transported  in  the  forward  and  the  backward  direc¬ 
tions  can  be  visualized  in  Figure  4,  which  shows  a  time  sequence  of  the  the  x-vorticity  contours  in  the 
v-2  plane  at  x  =  0.  Figure  4a  shows  the  instant  when  the  jet  first  impinges  on  the  ground.  Further  in 
time,  the  wall  jet  is  formed  in  both  the  forward  and  backward  directions.  However,  the  vorticity  on 
the  ring  behind  the  jet  is  much  lower  than  that  ahead  of  the  impingement  point.  The  size  of  the  vor¬ 
tex  behind  the  inclined  impinging  jet  is  also  smaller  than  that  ahead  of  the  jet  and  the  vortex  structure 
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Figure  3.  Three-dimensional  perspective  of  the  total  vorticity,)wj  in  the  inclined  impinging  jet 
for  the  unforced  case. 


c.  t  =  7.5 


Figure  4.  Time  sequence  of  the  propagation  of  the  starting  vortex  in  the  inclined  impinging  jet. 
The  .Y— vorticity  contours  in  the  y—z  plane  at  x  —  0.  Contour  interval  is  0.7. 
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appears  to  be  isolated.  These  figures  also  clearly  show  that  most  of  the  fluid  in  the  impinging  jet  is 
carried  in  the  forward  direction.  A  region  of  secondary  vorticity  (of  opposite  sign)  appears  ahead  of 
the  propagating  vortex  ring.  This  is  consistent  with  experimental  observation3  and  was  also  seen  in 
earlier  calculations.9  The  details  of  the  secondary  vorticity  could  not  be  determined  due  to  the  lack 
of  adequate  grid  resolution  in  the  wall  region. 

Figure  5  shows  the  corresponding  time  sequence  of  the  y-vorticity  contours  in  the  x-z  plane  at 
y  =  0.  This  is  a  view  in  the  plane  of  the  fountain.  The  shed  vortex  ring  propagates  toward  the  foun¬ 
tain  axis  (i.e.,  x  =  xT)  as  time  progresses.  If  the  computations  were  continued,  this  vortex  ring  in  the 
wall  jet  will  collide  with  the  vortex  ring  in  the  opposing  wall  t  and  create  an  upwash  in  the  fountain. 
This  upwash  would  then  spread  upward  as  the  vortex  ri  ..s  climb  up  the  fountain  and  eventually 
impinge  on  the  top  wall.  This  results  in  the  lift-off  effect  of  the  fountain  that  is  of  great  importance 
for  V/STOL  aircraft  operating  near  the  ground. 

Due  to  the  time  and  resource  constraints  of  the  Phase  I  study,  however,  the  simulation  was  not 
continued  until  a  steady  state  was  reached.  Instead,  the  solution  at  the  time  shown  in  Figure  3  was 
perturbed  by  adding  a  fluctuation  w ',  of  the  form  given  in  Equation  (11),  to  the  steady-state  jet  pro¬ 
file.  A  forcing  amplitude  of  0.2(7  was  employed  for  these  simulations.  Simulation  was  then  continued 
for  six  cycles  of  forcing  and  the  results  were  acquired  in  detail  and  analyzed  for  the  last  forcing  cycle. 

Figure  6a  shows  the  three-dimensional  perspective  of  the  total  vorticity  |w|  at  the  end  of  the  last 
forcing  cycle,  using  the  same  viewing  angle  and  vorticity  level  as  in  Figure  3.  Figures  6b  and  6c 
show  the  x-  and  y-vorticity  components  (i.e.,  |wt|  and  laj),  respectively,  in  a  similar  three-dimensional 
perspective.  The  bulge  in  the  impinging  jet  is  the  result  of  the  passage  of  the  shed  vortex  ring.  With  a 
higher  grid  resolution,  the  shed  vortex  ring  would  be  clearly  visible,  as  was  shown  in  earlier  calcula¬ 
tions.10  Due  to  the  forcing,  the  vortex  ring  is  shed  at  the  forcing  frequency.  The  impingement  of  this 
vortex  ring  on  the  ground  is  quite  similar  to  that  observed  for  the  starting  vortex  (Figure  3).  When  the 
vortex  ring  propagates  on  the  ground,  most  of  the  vorticity  is  carried  in  the  forward-moving  portion  of 
the  ring  while  the  vorticity  in  the  ring  behind  the  jet  is  much  weaker.  When  the  vortex  rings  reach  the 
fountain  axes,  they  collide  with  the  vortex  ring  from  the  adjacent  wall  jet  (due  to  the  periodic  boun¬ 
dary  conditions  at  x  jD  —  1.5).  After  the  collision,  the  vortex  rings  continue  to  propagate  up  into  the 
fountain.  The  lift-off  of  the  vortex  ring  into  the  fountain  can  be  seen  in  Figure  6a.  Due  to  the  com¬ 
bined  effect  of  the  lift-off  into  the  fountain  (in  the  x-z  plane)  and  the  forward  propagation  to  the  out¬ 
flow  (in  the  y-z  plane),  the  vortex  rings  starts  to  twist  and  eventually  breaks  down  into  smaller  struc¬ 
tures. 

The  vorticity  shed  by  the  impinging  jet  is  redistributed  due  to  the  motion  in  the  three-dimensional 
space,  as  can  be  observed  in  Figures  6b  and  6c.  For  example,  the  vorticity  in  the  outflow  direction  is 
predominantly  |w*|  (Figure  6b),  whereas  in  the  plane  of  the  fountain,  dominates.  The  region  of 
secondary  vorticity  on  the  ground  below  the  primary  vortex  ring  can  also  be  seen.  The  details  of  the 
small-scale  structures  in  the  flow  field  are  not  very  clear  in  these  figures,  primarily  due  to  the  low  grid 
resolution  of  the  calculations.  High-resolution  computations  of  normal  impinging  jets  show  more 
details  of  the  breakdown  of  the  large  structures.910 

Figure  7  shows  a  time  sequence  of  the  x-vorticity  contours  in  the  y-z  plane  at  x  =  0  for  the  entire 
last  forcing  cycle.  The  figures  are  shown  tt/2  apart  in  the  forcing  cycle.  In  comparing  the  first  and 
the  last  figure,  it  is  evident  that  the  sequence  of  vortex  shedding  in  the  jet  repeats  itself.  This  is 
characteristic  of  external  forcing,  which  suppresses  the  natural  instability  of  the  jet,  and  the  shed 
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a.  t  =  3.75  b.  t  =  5.63 


Figure  5.  Time  sequence  of  the  propagation  of  the  starting  vortex  in  the  inclined  impinging  jet. 
The  v— vorticity  contours  in  the  x—z  plane  at  v  =  0.  Contour  interval  is  0.5. 
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structures  are  organized  by  the  forcing  frequency.  The  phase  speed  of  the  vortex  ring  in  the  jet  shear 
layer  can  be  estimated  from  this  time  sequence  to  be  around  0.6 U,  which  agrees  with  the  experimental 
value  obtained  for  normal  impinging  jets.3  This  implies  that  the  obliqueness  of  the  impinging  jet  docs 
not  significantly  affect  the  convective  speed  of  the  structures  in  the  jet.  A  region  of  secondary  vorti- 
city  ahead  of  the  propagating  vortex  ring  is  again  observed.  It  appears  that  some  of  this  secondary 
vorticity  is  due  to  the  forward  transport  of  vorticity  from  the  rear  of  the  jet.  This  also  results  in  a 
reduction  of  the  vorticity  levels  behind  the  jet  (i.e.,  y/D  >  1).  The  details  of  the  secondary  vorticity 
arc  not  seen  clearly  in  this  simulation  due  to  the  coarseness  of  the  mesh.  In  the  earlier  study,®  it  was 
demonstrated  that  the  secondary  vorticity  appears  as  a  result  of  the  unsteady  separation  of  the  wall 
boundary  layer  caused  by  an  adverse  pressure  gradient.  This  behavior  has  also  been  determined  in 
experimental  studies.3  The  physics  of  this  phenomena  should  not  be  affected  by  the  inclination  of  the 
jet  and  thus  a  similar  process  of  unsteady  separation  is  expected  to  occur  for  the  inclined  jet  case.  To 
confirm  this,  however,  high-resolution  simulations  must  be  performed.  This  will  be  considered  in  the 
next  phase. 

The  separation  between  the  structures  in  the  wall  jet  give  a  wavelength  A  of  around  0.56D,  which 
for  a  forcing  frequency  /  of  0.47,  results  in  a  convective  speed  of  V  c  =  A/  =  0.26V.  This  is  again  simi¬ 
lar  to  that  observed  in  the  axisymmetric  forcing  study  of  normal  impinging  jets.®  This  speed  is  lower 
than  that  for  an  isolated  impinging  jet3  and  is  primarily  due  to  the  presence  of  the  fountain,  which 
retards  the  forward  motion  of  the  structures  in  the  outflow  direction. 

Figure  8  shows  the  corresponding  j-vorticity  contours  in  the  x-z  plane  at  y  =  0.  Due  to  the  sym¬ 
metry  conditions  at  .x  =  0,  only  half  of  the  full  domain  in  the  x-z  plane  is  shown.  This  view  shows  the 
motion  in  the  fountain  region.  Notice  that  the  vortex  ring  propagates  on  the  ground  and  is  then  lifted 
into  the  fountain.  The  vorticity  in  the  ring  decreases  as  it  moves  up  in  the  fountain,  primarily  due  to 
the  three-dimensional  spreading.  The  fountain  spreads  upward  and  impinges  on  the  undersurface  of 
the  top  plate  and  consequently  spreads  on  the  top  wall.  This  impingement  process  is  the  primary 
mechanism  for  the  lift-off  effect  of  the  fountain  on  a  V/STOL  aircraft.  Due  to  the  close  proximity  of 
the  jet  axis  to  the  fountain  axis  in  this  simulation,  some  portion  of  the  fluid  in  the  fountain  is 
entrained  into  the  impinging  jet.  This  would  result  in  a  modification  of  the  spreading  rate  of  the  imp¬ 
inging  jet. 

Figure  9  compares  the  velocity  vector  field  and  the  pressure  field  in  the  y-z  plane  at  the  end  of  the 
last  forcing  cycle.  The  inclined  impingement  of  the  jet  is  clearly  seen  (Figure  9a).  In  addition,  it  can 
be  observed  that  the  stagnation  point  on  the  ground  plane  does  not  occur  at  the  y  =  0  location  (which 
is  the  stagnation  point  for  a  normal  impinging  jet),  but  is  offset  ahead  by  a  certain  distance.  This  new 
location  of  the  stagnation  point  depends  upon  the  angle  of  impingement.  The  jet  physical  centerline  is 
inclined  at  6  to  the  ground  and  is  ahead  of  the  stagnation  streamline.  These  general  features  arc  in 
good  agreement  with  past  experimental  observations.14  Figure  9b  shows  the  velocity  vector  field  in 
the  plane  of  the  fountain  {x /D  =  1.5).  Unlike  the  symmetric  fan-like  behavior  of  the  upflow  in  the 
fountain  for  the  normal  impingement  case,  the  upwash  is  skewed  in  the  direction  of  inclination  of  the 
impinging  jet.  This  implies  that  the  flow  impingement  on  the  top  wall  will  also  be  skewed  and  thus 
will  probably  result  in  a  nonuniform  lift  force  on  the  aircraft.  Comparing  the  velocity  (Figure  9a)  and 
pressure  gradient  (Figure  9c)  vector  fields  with  the  corresponding  vorticity  contours  shown  earlier 
(Figure  7e)  indicates  that  there  is  an  adverse  pressure  gradient  ahead  of  the  primary  vortex  ring. 
Here,  the  direction  of  the  vector  indicates  the  direction  of  increasing  pressure.  There  is  also  a  low- 
pressure  region  associated  with  the  core  of  the  shed  vortex  ring.  These  low-pressure  regions  arc  seen 
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more  clearly  in  the  pressure  contours  in  the  y-z  plane  shown  in  Figure  9d.  The  pressure  contours  on 
the  ground  also  show  that  the  stagnation  pressure  is  offset  from  the  y  =  0  location.  This  distance  can 
be  determined  from  this  Figure  to  be  around  0.57 D.  The  pressure  contours  also  show  a  compression 
occurring  behind  the  peak  location.  This  behavior  was  also  observed  in  earlier  experiments5  and  in 
computations.16 

The  stagnation  streamline  can  be  approximately  located  by  using  the  stagnation  point  location  on 
the  ground  in  Figure  9a  or  9d  and  drawing  a  line  at  60'  to  the  ground.  Its  corresponding  location  in 
the  vorticity  plot  (Figure  7e)  can  then  be  determined.  It  can  then  be  shown  that  there  is  nonzero  vor- 
ticity  on  the  stagnation  streamline,  thereby  confirming  the  previous  predictions18  that  oblique 
impingement  results  in  a  stagnation  streamline  with  nonzero  vorticity. 

Figures  10a  and  10b  show  the  velocity  vector  field  in  the  x-z  plane  at  y  =  0  and  y  /D  =  1.5. 
respectively.  The  comparison  of  the  magnitude  of  the  velocity  vectors  in  the  impinging  jet  (Figure 
10a)  and  in  the  fountain  (Figure  10b)  indicates  that  the  fountain  is  relatively  weaker.  This  is  due  to 
obliqueness,  which  results  in  a  significant  velocity  component  in  the  positive  y-direction  and  thus 
changes  the  direction  of  fluid  motion  in  the  fountain  region.  Some  of  the  fluid  impinging  on  the  top 
wall  spreads  and  is  entrained  into  the  impinging  jet.  This  spreading  and  the  entrainment  process  arc 
stronger  away  from  the  jet  centerline,  as  seen  in  Figure  10b.  Recirculation  of  the  fluid  occurs  at  this 
location  and  is  quite  strong  due  to  the  short  jet-to-ground  height.  If  the  height  H  is  increased,  this 
recirculation  would  be  weakened  and  the  effects  of  entrainment  of  the  fluid  into  the  impinging  jet 
would  also  be  reduced. 

Figure  1 1  shows  the  pressure  contours  on  the  top  plate  and  the  ground  plane.  The  offset  between 
the  center  of  the  jet  at  its  exit  plane  on  the  top  wall  (which  is  a  low-pressure  region)  and  the  stagna¬ 
tion  point  on  the  ground  is  clearly  observable.  As  seen  earlier,  the  pressure  field  on  the  ground  is 
skewed  in  the  direction  of  the  forward  motion  with  some  compression  occurring  behind  the  peak 
location.  It  is  possible  to  estimate  the  net  lift-off  on  the  top  plate  by  integrating  the  pressure  on  the 
plate.  From  this  figure,  the  offset  of  the  stagnation  point  for  the  inclined  jet  from  its  location  for  a 
normal  impingement  (at  y  =  0)  can  be  estimated.  A  value  of  around  0.57 D  is  again  estimated.  This 
computed  offset  compares  quite  well  with  the  experimental  data  of  Donaldson  and  Snedeker5  who 
obtained  a  value  of  around  0.588D  for  H  jD  =  1.96. 

Figure  12  shows  typical  velocity  profiles  in  the  flow  field.  Figure  12a  shows  the  v-velocity  profiles 
in  the  y-z  plane  at  x  =  0.  The  profiles  one  half-cycle  apart  in  the  last  forcing  cycle  are  shown  in  this 
figure.  It  is  clear  that  in  the  forward  direction  (i.e.,  the  negative  y-direction)  the  wall  boundary  layer 
grows  toward  the  outflow.  However,  in  the  backward  direction  (positive  y-direction),  the  wall  jet  is 
quite  weak,  and  very  little  outflow  is  occurring  near  the  outflow  boundary.  This  also  indicates  that 
most  of  the  fluid  is  being  transported  in  the  forward  direction.  The  effect  of  forcing  is  observable 
only  near  the  wall  region. 

Figure  12b  shows  the  corresponding  u-velocity  profile  in  the  x-z  plane  at  y  =  0.  Due  to  the  close 
proximity  of  the  fountain  axis,  the  effect  of  forcing  is  quite  noticeable  in  these  velocity  profiles.  The 
growth  of  the  wall  jet  is  again  observable.  On  the  top  wall,  there  is  a  reverse  wall  jet  which  is  caused 
by  the  impingement  of  the  fountain  on  the  top  wall.  The  reverse  wall  jet  also  causes  the  fluid  in  it  to 
be  entrained  into  the  impinging  jet  and  thus  modifies  the  spreading  rate  of  the  impinging  jet. 
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a.  v  =  0  plane. 


Figure  10.  The  velocity  vector  fields  in  the  v—z  plane  at  t= 
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.  Top  wall. 


b.  Ground  plane. 


Figure  11.  Pressure  contours  in 
the  x-y  plane  at  t=24.375. 


c.  vv— velocity  profiles  in  the  y—  z  plane  at  x  =  0. 
Figure  12.  The  velocity  profiles  in  the  flow  field  (cont.). 


Figure  13.  The  pressure  variation  on  the  ground  plane. 


Figure  12c  shows  the  w-velocity  component  in  the  y-z  plane  at  x  =  0.  This  figure  show-s  the  verti¬ 
cal  velocity  profiles  in  the  impinging  jet.  The  spreading  of  the  impinging  jet  and  the  oblique  impinge¬ 
ment  of  the  jet  can  be  visualized  from  this  figure.  The  effect  of  forcing  on  the  profiles  near  the 
ground  is  also  seen.  Notice  that  there  is  some  positive  >v-velocity  close  to  the  wall.  This  implies  an 
upflow  and  is  indicative  of  lift-off  from  the  ground  associated  with  the  separation  and  the  formation 
of  secondary  vorticity  ahead  of  the  propagating  vortex  rings  in  the  wall  jet.  A  similar  behavior  was 
observed  in  the  earlier  normal  impingement  study.9 

In  nrinciple,  other  data  such  as  the  jet  centerline  velocity  decay  could  be  reduced  from  the  com¬ 
puted  results.  However,  since  the  computational  grid  is  Cartesian  and  the  jet  centerline  is  inclined,  a 
complex  three-dimensional  interpolation  of  the  centerline  velocity  field  from  the  computed  velocities 
on  the  computational  nodes  is  required  to  estimate  the  centerline  velocity.  Due  to  the  limited 
resources  of  the  Phase  I  study,  the  development  of  the  interpolation  package  has  been  deferred  to  the 
next  phase. 

Finally,  Figure  13  shows  the  variation  of  the  pressure  on  the  ground  plane  in  the  y-z  plane  at 
,v  =  0.  The  pressure  at  the  outflow  boundary  at  y  /D  =  3.0  for  the  normal  impingement  case  was 
assumed  to  represent  the  ambient  pressure  (Poo)-  The  pressure  is  normalized  by  the  pressure  for  the 
normal  impingement  case  (to  be  discussed  in  Section  5.2.1)  by  the  relation 


P 


P  -  Poo 
( Po)qO  P  oo 


(12) 


where  (pJoo  is  the  stagnation  pressure  for  the  normal  impingement  case.  This  normalization  was  car¬ 
ried  out  to  be  consistent  with  the  experimental  normalization.5  The  pressure  variation  for  the  normal 
impingement  is  shown  in  this  figure  along  with  the  pressure  variation  for  the  inclined  jet  at  two 
instances  in  the  forcing  cycle.  Only  qualitative  comparison  between  the  normal  and  the  oblique 
impingement  is  intended  here  since  the  simulation  with  normal  impingement  (described  in  the  next 
section)  was  not  carried  out  to  steady  state.  Thus,  the  higher  value  of  the  peak  pressure  for  the  nor¬ 
mal  impingement  case  is  probably  due  to  the  transient  effects  in  that  simulation.  The  shift  of  the  stag¬ 
nation  pressure  due  to  oblique  impingement  can  be  observed  in  this  figure.  The  effect  of  forcing 
causes  the  peak  pressure  to  vary  with  the  forcing  cycle.  However,  the  location  is  not  changed  signifi¬ 
cantly.  Note  also  that,  whereas  the  pressure  distribution  is  symmetric  for  the  normal  impingement 
case,  oblique  impingement  results  in  significant  asymmetry.  There  is  an  unsteady  variation  of  pres¬ 
sure  away  from  the  stagnation  point,  with  the  forcing  causing  the  favorable  pressure  gradient  (at 
t  =  22.97)  to  become  an  adverse  pressure  gradient  (at  t  =  24.375). 


5.2  The  Normal  Impinging  Jet  without  Crossflow 

Since  the  normal  impinging  jet  was  also  simulated,  a  brief  discussion  of  the  results  for  this  case  is 
also  presented  here  for  completeness.  Furthermore,  the  effect  of  oblique  impingement  (Section  5.1) 
and  the  effect  of  crossflow  (to  be  discussed  in  Sections  5.3  and  5.4)  can  be  understood  in  greater  detail 
through  comparison  with  the  reference  case  of  a  normal  impinging  jet.  More  detailed  analyses  of  the 
normal  impinging  jets  have  been  described  earlier.910  In  the  present  study,  the  simulation  was  carried 
out  only  for  a  short  duration  and  the  steady-state  solution  was  not  reached.  Thus,  only  a  qualitative 
interpretation  of  the  computed  results  is  presented  here. 
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The  simulation  for  the  normal  impingement  case  was  also  begun  with  no  flow  and  with  the  jet 
velocity  profile  given  in  Equation  (10)  imposed  at  the  jet  exit  at  t  =  0.  Figure  14  shows  a  three- 
dimensional  perspective  of  the  total  vorticity  M  for  the  normal  impingement  case  at  the  end  of  the 
present  simulation.  At  this  time,  the  starting  vortex  has  just  reached  the  fountain  axis  and  has  begun 
to  lift  off.  The  flow  field  is  symmetric  in  both  the  x-z  and  y-z  planes,  as  expected,  and  the  secondary 
vorticity  ahead  of  the  primary  vortex  ring  can  be  seen.  Although  the  perspective  viewing  angle  is  dif¬ 
ferent  from  that  shown  for  the  oblique  impinging  case  (Figure  3),  comparison  between  the  figures 
shows  the  characteristic  changes  due  to  oblique  impingement. 

Figures  15a  and  15b  show  the  x-vorticity  in  the  y-z  plane  at  x  =  0  and  the  y-vorticity  in  the  x-r 
plane  at  y  =  0,  respectively,  at  the  same  point  in  time  shown  in  Figure  14.  The  symmetric  propaga¬ 
tion  of  the  vortex  ring  in  the  wall  jet  (Figure  15a)  and  the  lift-off  into  the  fountain  in  the  x-z  plane 
(Figure  15b)  can  be  seen.  Notice  that,  due  to  the  proximity  of  the  fountain  axes,  the  vortex  ring  rises 
up  in  the  upwash  while  in  the  outflow  direction  it  is  still  propagating  on  the  ground.  This  eventually 
results  in  the  breakdown  of  the  ring  into  smaller  structures. 

Figure  16  shows  representative  velocity  vector  fields  in  the  flow  field.  Figures  16a  and  16b  show 
the  velocity  fields  in  the  y-z  plane  at  x  =  0  and  at  x  /D  =  1.5,  respectively.  The  vortical  motion  in  the 
wall  jet  causes  the  flow  to  rotate  around  the  vortex  ring  and  entrain  the  potential  field  into  the  wall  jet 
(Figure  16a).  The  stagnation  point  of  the  impinging  jet  is  located  at  y  =  0,  as  expected.  Comparing 
Figures  15a  and  16a,  it  can  be  seen  that  the  vortex  2nd  the  velocity  field  are  not  at  exactly  the  same 
location.  As  shown  previously,9  it  is  incorrect  to  infer  the  vortical  motion  (which  is  the  result  of  the 
velocity  gradients)  from  the  velocity  field  alone.  This  is  an  important  factor  in  interpreting  the  flow 
field  associated  with  vortical  motion. 

Figure  16b  shows  the  symmetric  fan  in  the  plane  of  the  fountain.  The  upwash  has  not  yet  reached 
the  top  wall  at  this  stage  of  the  computations.  This  figure  can  be  compared  with  the  corresponding 
case  for  oblique  impingement  (Figure  9b)  to  identify  the  differences  in  the  upwash  due  to  oblique 
impingement.  Figure  16c  shows  the  velocity  vector  field  in  the  x-z  plane  at  y  =  0.  The  slowing  of  the 
wall  jet  as  the  fountain  axis  is  reached  can  be  observed.  There  is  another  stagnation  point  at  the  loca¬ 
tion  where  the  two  wall  jets  collide  that  is  slightly  above  the  ground,  as  can  also  be  seen  in  Figure 
16b.  The  formation  of  the  upwash  is  also  apparent  in  this  figure.  Comparison  with  the  corresponding 
plot  for  the  inclined  jet  (Figure  10a)  shows  that  there  is  more  recirculation  occurring  in  the  normal 
impingement  case.  Thus,  for  this  H  /D  case,  the  entrainment  of  the  fluid  into  the  impinging  jet  due  to 
this  recirculation  is  much  stronger  for  the  normal  impingement  case. 

Finally,  Figures  17a  and  17b  show  the  pressure  contours  in  the  x-y  plane  at  z  /H  =  1  and  z  /H  -  0. 
respectively.  On  the  top  plate,  the  pressure  field  is  nearly  uniform  in  the  core  of  the  jet  and  the  pres¬ 
sure  increase  in  the  region  where  the  fountain  impinges  on  the  top  wall  can  be  seen.  On  the  bottom 
wall,  the  stagnation  point  is  on  the  jet  centerline,  as  mentioned  earlier.  Also,  the  pressure  field  is  sym¬ 
metric,  unlike  the  oblique  impingement  case  in  which  the  pressure  field  on  the  ground  was  skewed, 
with  compression  occurring  behind  the  stagnation  point  (Figure  1  lb). 

The  results  presented  above  show  the  characteristic  flow  field  associated  with  the  normal  impinge¬ 
ment  case.  Simulations  of  excited  normal  impinging  jets  carried  out  earlier910  also  showed  complex 
vortical  motion.  In  the  following  sections,  we  discuss  the  effect  of  crossflow  on  the  normal  impinging 
jet. 
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c.  x  plane  at  y  =  0. 


Figure  16.  The  velocity  vector  Helds  for  the  normal 
impinging  jet  at  t  =  9.4. 


a.  Top  wall. 


b.  Cround  plane. 


Figure  17.  The  pressure  contours  in  the 
x  — v  plane  at  t  =  9.4. 
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5.3  The  Normal  Impinging  Jet  in  Crossflow 

The  effect  of  crossflow  on  the  dynamics  of  fluid  motion  in  the  impinging  jet  was  studied  by  impos¬ 
ing  a  crossflow  velocity  at  the  inflow  boundary  at  y/D  =  -3.0.  Two  types  of  simulations  were  per¬ 
formed.  In  the  first  simulation,  the  solution  obtained  at  the  end  of  the  simulation  with  the  normal 
impinging  jet  (Section  5.2)  was  used  to  initialize  the  flow  field  and  then  the  crossflow  was  introduced. 
The  second  simulation  was  started  with  no  flow  in  the  computational  domain  and  with  the  crossflow 
and  the  jet  profile  given  in  Equation  (10)  imposed  simultaneously  at  t  =  0.  In  terms  of  the  flight  con¬ 
figuration  of  a  V/STOL  aircraft  in  ground  effect,  these  two  scenarios  correspond  approximately  to  the 
cases  in  which  (a)  the  aircraft  executes  a  forward  motion  after  hovering  above  the  ground,  and  (b)  the 
aircraft  moves  forward  as  soon  as  it  lifts  off  the  ground.  Although  the  simulations  were  not  carried 
out  for  a  sufficient  length  of  time  to  reach  a  quasi-stationary  state,  the  results  shown  below  indicate 
that  the  flow  field  associated  with  these  two  flight  situations  can  be  quite  different.  A  constant 
crossflow  velocity,  vQ  =  Vc/U  =  0.4,  where  Vc  is  the  crossflow  velocity,  was  used  in  all  the  simulations. 
The  effect  of  varying  the  crossflow  velocity  will  be  studied  in  more  detail  in  the  next  phase.  Due  to 
the  computer  and  resource  restrictions  of  the  Phase  I  study,  the  effect  of  forcing  was  studied  only  for 
the  latter  case  and  only  a  few  forcing  cycles  were  performed. 

Note  that  the  presence  of  the  fountain  modifies  the  flow  field  significantly,  a  phenomenon  that  is 
absent  in  most  of  the  earlier  experimental  and  numerical  studies  of  jets  in  crossflow. 

5.3.1  Forward  Motion  after  Hovering  (Case  A) 

To  simulate  this  configuration,  the  crossflow  was  imposed  at  the  inflow  boundary  at  y/D  =  -3.0  on 
the  solution  obtained  for  the  normal  impingement  case  described  in  Section  5.2.  The  crossflow  was 
increased  from  v0  =  0  to  v0  =  0.4  within  25  time  steps,  thereby  simulating  a  rapid  initial  acceleration 
from  stationary  hover  to  forward  motion  at  a  constant  velocity  v0.  The  simulation  was  continued 
without  any  external  excitation  and  the  data  were  obtained  in  detail  at  the  end  of  the  simulation.  Due 
to  time  and  resource  constraints,  the  computation  was  not  carried  out  for  a  sufficient  length  of  time  to 
reach  a  steady  state.  However,  analysis  of  the  data  indicated  that  the  solution  was  converging  and 
w'ould  eventually  reach  a  steady  state.  Some  representative  results  are  presented  and  analyzed  here. 

Figure  18a  shows  the  three-dimensional  perspective  of  the  total  vorticity  |w|  in  the  flow  field  at  the 
end  of  this  simulation.  The  jet  had  already  impinged  on  the  ground  and  the  lift-off  of  the  vortex  ring 
into  the  fountain  had  begun  (Figure  14)  when  the  crossflow  was  initiated.  The  crossflow  quickly 
forces  back  the  forward  part  of  the  vortex  ring  on  the  ground  (cf.  Section  5.2,  Figure  15a)  and  wraps 
the  structure  around  the  jet,  since  the  front  part  of  the  impinging  jet  acts  as  an  obstacle  to  the  cross 
stream  and  thus  the  crossflow  must  go  around  it.  The  rear  portion  of  the  jet  is  also  spread  in  the  posi¬ 
tive  ^'-direction  due  to  the  crossflow. 

To  visualize  the  flow  field  in  more  detail,  Figures  18b  and  18c  show  the  |«x|  and  la^l  components, 
respectively,  in  the  same  three-dimensional  perspective.  The  |wt|  component  dominates  the  vorticity 
in  the  ring  ahead  of  the  jet  and  in  the  ring  behind  it  on  the  ground  plane,  whereas  the  vorticity  in  the 
fountain  regions  is  dominated  by  the  \w^  component.  It  has  been  shown  in  earlier  experimental  stu¬ 
dies  (e.g.,  Andreopoulos  and  Rodi4)  that  the  dynamics  of  the  vorticity  redistribution  is  quite  complex 
for  jets  in  crossflow.  Further  interpretation  of  the  vorticity  redistribution  mechanism  is  given  in  some 
detail  in  the  next  section. 
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Figures  19a  and  19b  show  the  .x-vorticity  in  the  y-z  plane  at  x  =  0  and  the  v-vorticity  in  the  v-r 
plane  at  y  =  0,  respectively.  The  forward-moving  part  of  the  original  vortical  ring  (Figure  15a)  has 
been  pushed  back  due  to  the  crossflow  and  is  now  embedded  near  the  jet  impingement  region.  The 
jet  itself  is  now  inclined  in  the  positive  y-direction  and  significant  secondary  vorticity  (dotted  con¬ 
tours)  is  now  present  below  the  primary  vortical  region  in  the  wall  jet  for  y  /D  >  0.  Part  of  the  secon¬ 
dary  vorticity  in  this  region  is  now  due  to  the  vorticity  in  the  ring  ahead  of  the  jet  being  forced  back 
by  the  crossflow.  In  the  plane  of  the  fountain  (Figure  19b),  the  vortical  structure  seen  in  the  fountain 
axis  for  the  normal  impinging  case  (cf.  Figure  15b)  is  now  located  away  from  the  fountain  axis.  This 
is  another  effect  of  the  crossflow  on  the  fountain  which  results  in  adding  a  v-component  of  the  velo¬ 
city  to  the  primarily  w-component  of  velocity  in  the  fountain  upwash.  The  structures  in  the  plane 
v/D  =  1.5  are  more  complex,  as  seen  in  Figure  19c.  In  this  figure,  note  that  the  vorticity  near  the 
ground  is  associated  with  the  diffused  jet  shear  layer  which  has  spread  in  the  positive  y-direction  due 
to  crossflow.  The  vortical  structures  near  the  top  wall  are  part  of  the  cross-sectional  view  of  the  vor¬ 
tex  ring  in  the  fountain.  Figures  18  and  19  should  be  used  together  to  visualize  the  complex  flow 
field. 

Figure  20  shows  some  representative  velocity  vector  fields  and  the  pressure  field.  Figure  20a 
shows  the  velocity  field  in  the  y-z  plane  at  x  *=  0  corresponding  to  the  vorticity  field  shown  In 
Figure  19a.  The  crossflow  at  v/D  =  -3.0  is  clearly  seen  and  the  formation  of  the  recirculation  ahead 
of  the  jet  is  also  observable.  This  recirculation  is  associated  with  the  vortex  seen  in  the  vorticity  field 
(Figure  19a).  Due  to  the  addition  of  the  crossflow  component,  the  outflow  velocity  in  the  wall  jet 
region  is  increased.  Notice  also  that  the  stagnation  point  is  moved  downstream  of  the  y  =  0  location 
where  it  originally  was  before  the  crossflow  was  started.  The  symmetric  fan-like  upwash  for  the  nor¬ 
mal  impingement  case  (Figure  16b)  at  the  plane  of  the  fountain  at  x  /£>  =  1.5  is  no  longer  present  with 
crossflow,  as  can  be  seen  in  Figure  20b.  The  upflow  is  now  skewed  in  the  positive  y-direction  due  to 
the  addition  of  the  crossflow  velocity  to  the  upwash  velocity,  causing  the  flow  to  turn  downstream. 
This  would  result  in  a  much  weaker  lift-off  effect  (if  any)  on  the  top  wall. 

The  pressure  contours  in  the  y-z  plane  shown  in  Figure  20c  illustrate  the  skewness  in  the  positive 
y-direction  due  to  crossflow.  There  is  a  low-pressure  region  associated  with  the  vortex  ahead  of  the 
jet.  The  stagnation  point  is  no  longer  on  the  jet  center  line  and  is  moved  by  about  0.27D  to  the  right 
of  the  y  =  0  location. 

Figures  20d  and  20e  show  the  velocity  vector  field  in  the  .x-c  plane  at  y  =  0  and  y /D  =  1.5. 
respectively.  The  recirculation  observed  for  the  normal  impingement  case  (Figure  16c)  is  no  longer 
present  in  the  center  plane,  and  the  upwash  velocity  in  the  fountain  is  considerably  weaker.  This  is 
not  surprising  since,  as  mentioned  earlier,  the  crossflow  in  the  y-direction  enforces  a  net  velocity  in 
that  direction  which  is  out  of  the  frame  of  this  figure.  Halfway  to  the  outflow  at  v/D  =  1.5,  we  see 
significant  recirculation  due  to  the  complex  vortex  motion  observed  in  Figure  19c.  The  upwash  in  the 
fountain  region  is  quite  strong  in  this  plane. 

Figure  21  shows  some  representative  velocity  profiles  in  the  flow  field.  For  comparison,  the  velo¬ 
city  profiles  for  the  normal  impingement  case  are  also  shown  in  these  figures.  Figure  21a  shows  the 
v-componcnt  profiles  in  the  y-z  plane  at  x  =  0.  For  the  normal  impingement  case,  the  velocity  profiles 
arc  symmetric  about  the  y  =  0  location,  with  the  wall  jet  growing  uniformly  toward  both  outflow 
boundaries.  With  crossflow,  the  profiles  at  all  the  locations  show  motion  primarily  in  the  positive  y- 
dircction.  The  wall  boundary  layer  due  to  the  crossflow  can  be  seen  in  this  figuie.  The  incoming 
crossflow  velocity  profile  in  the  potential  flow  region  is  uniform  across  the  entire  H /D  at  v/D  -  -2.33 
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a.  x—  vorticity  in  the  y—z  plane  at  x  =  0 


i  i  i 

. . . .  ¥ 1  i  V  V*f? 1  ii 

i  i  i  mm  r  i  IV*',  i  [ 

i  . . .  i  r"i 

i  I  IIIIIIIII  yi  'l 
I  I  mum/,1 Bill 

. mi'  S'  i 

i  i  mini //H  i 
i  i  innn//  i  '  M 

/ 1 1 1  n  / 1©  i  i 

/ 1  iii  v  //A7'/ 

iiin  \ ,7,71, 

I  I'M 

'I'll  I  II  I  I 

It'll  si  l  I 

'  1 :  V7/  / : 

I  v, 

•  -.'ii 

I  V  //  - 
1 

\  /  / 


"/  'i  ]'; 

/V7//<'/7/(/  i  c  /  J  1  1 

//'/I*  \  N  V  - - "  / 

^  /  /  /  /  ,  i  i  ,  \  ^ _ ✓  ;  i 

<v  /  /  i  i  \  \ 

/  '  I  v  I  I  \  --  / 

'  I  I  >  V  -1. 60-  i 


- .  \  \  \ '  1 1 1 , 1 

rr  / ;  i  'jin 

L  /  'll  III 

I  !  Il|l  1 

^  )  i  nuJ’ 

--  - /  nS1 

--'-1.60-  nr?1 

/'I1! 

in  1 


;  i  / 

i  i  i 

I,/  s 

i  i  § 

i  i  ' 

I  v  I  <M 

\  \  \  ^ 
\  \  V 

\  s 


/  /  it  ' 

1,11  I  I 

/  /  I  / 


\  \ g  / '  ml' 

' . \  i  i;'///\  |J[!! 

1  '  1 1 1  \  &i 

- - _  _ i  "i'J.'i 


TR-470/02-89 


34 


b.  Velocity  vectors  in  the  y—  z  plane  at  x  /D  =  1.5 


c.  Pressure  contours  in  the  y—z  plane  at  x  =  0. 


Figure  20.  The  velocity  vector  fields  and  pressure  contours  for  normal  impinging  jet 
crossflow  (case  A)  at  t  =  13.1. 
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c.  Velocity  vectors  in  the  .v-r  plane  at  v  /D  =  1.5. 


Figure  20.  The  velocity  vector  fields  and  pressure  contours  for  normal  impinging  jet  with 
crossflow  (case  A)  at  t  =  13.1  (Cont.) 


TR-470/02-89 


36 


but  becomes  nonuniform  with  increasing  y  /D.  This  is  the  result  of  the  reversal  of  the  fluid  in  the  wall 
jet  in  the  region  y  /D  <  0  to  the  positive  y-direction.  The  wall  boundary  layer  grows  larger  for  the 
case  with  crossflow  for  y  /D  >  1.5  since  now  the  crossflow  velocity  component  is  added  to  the  velo¬ 
city  in  the  wall  boundary  layer.  All  these  features  occurred  as  expected  and  clearly  show  the  effect  of 
crossflow. 

Figure  21b  shows  the  »-velocity  component  in  the  x-z  plane  at  y  =  0.  In  this  plane,  the  profiles 
are  quite  different.  The  normal  impingement  velocity  profiles  are  now  more  developed  than  in  the 
case  with  crossflow,  with  significant  reverse  flow  near  the  top  wall  due  to  the  upwash  spreading  effect. 
With  crossflow,  however,  most  of  the  fluid  is  forced  to  move  to  the  outflow  (y-z  plane)  direction  and 
thus  the  motion  toward  the  fountain  region  in  the  x-z  plane  is  reduced.  This  then  results  in  reduced 
n-velocities  as  seen  in  this  figure.  Consequently,  the  upward  motion  in  the  fountain  is  weakened  and 
the  reverse  wall  jet  on  the  top  wall  is  also  reduced. 

Finally,  Figure  21c  shows  the  w-velocity  in  the  y-z  plane  at  x  =  0.  The  velocity  profiles  in  the  jet 
for  the  normal  impingement  case  is  again  symmetric  and  the  spreading  of  the  jet  is  clearly  seen.  Near 
the  initial  jet  exit  region  (z  /H  >  0.75),  the  effect  of  crossflow  is  insignificant  since  the  jet  peak  velo¬ 
city  is  2.5  times  the  crossflow  velocity.  As  the  jet  approaches  the  ground  plane,  however,  the  peak 
velocity  in  the  jet  potential  core  decreases  as  it  spreads  and  thus,  closer  to  the  ground,  the  crossflow 
velocity  is  of  the  same  order  as  the  w-component.  This  causes  the  jet  to  bend  toward  the  positive  y- 
direction,  and  the  spreading  character  of  the  jet  is  modified  by  the  crossflow.  The  crossflow  also 
forces  the  jet  to  impinge  a  short  distance  downstream  of  the  normal  impingement  location. 

In  this  section,  we  have  described  some  of  the  effects  of  crossflow  on  the  fluid  motion  of  a  normal 
impinging  jet.  The  characteristic  changes  in  the  features  of  vortical  and  fluid  motion  due  to  crossflow 
were  highlighted  through  comparison  with  the  reference  case  of  normal  impingement  without 
crossflow.  Although  the  simulations  were  not  carried  out  for  a  sufficient  length  of  time  to  reach  a 
quasi-stationary  state,  the  general  features  of  the  initial  modification  of  the  flow  field  are  in  good  qual¬ 
itative  agreement  with  experimental  observations. 

5.3.2  Forward  Motion  without  Hovering  (Case  B) 

For  this  flight  configuration,  the  crossflow  was  begun  at  t  =  0,  which  is  when  the  jet  comes  out  of 
the  jet  exit  plane  at  z  /H  =  1.  The  crossflow  was  again  ramped  to  the  steady  value  as  before.  Since 
the  crossflow  velocity  is  relatively  weaker  than  the  jet  velocity  and  since  the  jet-to-ground  distance  is 
quite  small  for  this  simulation,  the  jet  still  impinges  on  the  ground.  However,  the  flow  field  for  this 
case  develops  in  quite  a  different  manner  than  in  the  simulation  discussed  in  Section  5.3.1.  These 
differences  are  highlighted  below. 

Figure  22a  shows  the  three-dimensional  perspective  of  the  total  vorticity  in  the  jet  at  the  time  the 
jet  has  reached  the  ground.  There  is  no  clear  indication  of  the  starting  vortex  observed  in  the  earlier 
simulation,  and  the  jet  shear  layer  is  now  wrapped  behind  the  jet.  This  can  be  seen  more  clearly  in 
Figures  22b  and  22c,  which  show  the  x-vorticity  in  the  y-z  plane  at  x  =  0  and  the  y-vorticity  in  the 
x-z  plane  at  y  =  0,  respectively.  Due  to  the  crossflow,  the  jet  shear  layer  is  inclined  in  the  front  and 
more  spread  out  in  the  back  (Figure  22b).  The  vortex  seen  ahead  of  the  jet  in  Figure  19a  is  no  longer 
present  since  now  the  vortex  ring  does  not  exist  on  the  ground.  There  is  also  no  indication  of  the 
fountain,  which  was  expected  since  the  wall  jet  spreading  on  the  ground  plane  has  not  yet  occurred. 
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■Without  Crossflow 


- With  Crossflow 


0  .5 


-S  -2 


c.  n  —  velocity  profiles  in  the  v— :  plane  at  .v  =  0. 

Figure  21.  The  comparison  between  the  velocity  profiles  in  the  flowfield  for  the  normal 
impinging  jet  at  t  =  9.4  and  the  normal  impinging  jet  w  ith  crossflow  at  t  =  13.1.  (Cont.) 
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a.  3D  total  vorticity,|u/|at  a  level  of  1.3 

Figure  22.  The  vorticity  fields  for  the  normal  impinging  jet  with  crossflow  (case  B)  at  t  =  3.8. 
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b.  x— vorticity  in  the  y—z  plane  at  x  =  0.  Contour  interval  is  0.5 
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c.  y  — vorticity  in  the  x— z  plane  at  y  =  0.  Contour  interval  is  0.5. 

Figure  22.  The  vorticity  Helds  for  the  normal  impinging  jet  with  crossflow  (case  B)  at  t  =  3.8  (cont.) 


a.  x  =  0  plane 


b.  x  /£)  =  1.5  plane. 


Figure  23.  The  velocity  vector  fields  in  the  y—z  plane  for  normal  impinging  jet  with  crossflo" 
(case  B)  at  t  =  3.8. 
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The  bending  of  the  jet  due  to  crossflow  seen  in  Figure  22b  is  also  clearly  seen  in  the  velocity  vec¬ 
tor  field  in  the  same  plane  shown  in  Figure  23a.  Although  not  very  clear  in  this  figure,  there  is  a  small 
recirculation  region  ahead  of  the  jet  near  the  impingement  point.  The  velocity  vector  field  in  the 
plane  of  the  fo'mtaJn  at  x /D  =  1.5  is  shown  in  Figure  23b.  This  fig'"e  shows  that  the  fluid  flov.  in 
this  plane  is  primarily  in  the  outflow  direction  with  no  significant  upwash. 

The  computations  were  then  continued  with  external  excitation  at  the  same  Strouhal  number  of 
0.47  as  before.  A  forcing  amplitude  of  0.1  U  was  used  and  the  computations  were  continued  for 
approximately  four  forcing  cycles.  The  data  for  the  last  forcing  cycle  were  then  analyzed  in  detail. 

Figure  24a  shows  the  three-dimensional  perspective  of  the  vorticity  in  the  jet  at  the  end  of  the 
forced  simulation.  Comparison  with  Figure  18  from  the  earlier  study  of  the  crossflow  immediately 
reveals  significant  differences.  For  example,  the  fountain  does  not  form  at  all.  The  motion  toward  the 
fountain  axis  is  completely  inhibited  by  the  crossflow  and  most  of  the  fluid  is  swept  out  to  the  outflow 
direction.  Therefore,  in  this  case  there  will  be  no  lift-off  effect  of  the  fountain.  There  is  a  clear  indi¬ 
cation  of  the  counterrotating  vortical  structure  formed  by  the  jet  due  to  cross  flow.  This  is  quite  simi¬ 
lar  to  the  kidney-shaped  structures  observed  for  free  jets  in  crossflow,4’22  41  which  have  been  studied 
in  great  detail  in  the  past.  As  discussed  below,  however,  the  close  proximity  of  the  ground  plane 
modifies  the  shape.  The  wrapping  of  the  vortical  structure  around  the  jet  is  also  seen  in  this  figure. 
This  feature  is  similar  to  that  observed  in  the  earlier  simulation  (see  Figure  18). 

To  visualize  the  streamwise  (IwJ)  vortices  formed  due  to  crossflow.  Figure  24b  shows  the  \uiy\  com¬ 
ponent  in  the  same  perspective  view.  Comparing  this  figure  to  Figure  24a,  which  shows  the  total  vor¬ 
ticity,  it  becomes  quite  apparent  that  the  \u>^  component  of  the  vorticity  dominates  the  entire  flow 
field.  This  is  in  contrast  to  the  earlier  simulation,  where  the  presence  of  the  fountain  results  in  redis¬ 
tribution  of  the  total  vorticity  into  both  the  x-  and  .v-components  (Figures  18b  and  18c). 

To  illustrate  the  vortex  motion  in  the  jet.  Figures  25a  and  25b  show  the  x-vorticities  in  the  y-r 
plane  at  x  =  0  half  a  forcing  cycle  apart.  The  impinging  jet  is  swept  back,  similar  to  that  observed  in 
the  earlier  simulation  (Figure  19a).  Vortex  shedding  again  occurs  at  the  forcing  frequency;  however, 
the  motion  in  the  shear  layer  facing  the  crossflow  is  not  very  clear.  The  effect  of  the  crossflow  is  to 
compress  the  shear  layer  in  the  direction  of  the  crossflow,  and  some  of  the  vorticity  in  the  jet  shear 
layer  facing  the  crossflow  is  carried  into  the  wall  jet.  The  vortical  ring  that  is  shed  due  to  forcing  is 
highly  skewed  to  the  outflow  direction  and  is  elongated  in  the  wall  jet  region.  Notice  the  apparent 
termination  of  the  wall  jet  around  y  /D  =  2.  This  is  not  actually  occurring  and  is  due  to  the  fact  that 
further  downstream  the  wt  component  is  much  weaker  than  the  contour  interval  used  for  this  figure. 
Furthermore,  as  shown  earlier  (Figure  24b),  most  of  the  vorticity  is  contained  in  the  wv  component 
and  thus  rapid  depletion  of  the  w,  component  is  not  surprising.  In  fact,  it  appears  that  the  three- 
dimensional  redistribution  of  the  total  vorticity  in  the  flow  field  is  quite  complex,  as  seen  in  experi¬ 
ments,4  and  more  detailed  simulations  are  required  to  understand  the  detailed  dynamics  of  the  vorti¬ 
cal  motion  in  this  type  of  flow  field. 

The  y-vorticity  in  the  x-z  plane  at  y  =  0  and  y  /D  =  1.5  is  shown  in  Figures  25c  and  25d,  respec¬ 
tively.  Notice  the  complete  absence  of  a  fountain  in  the  y  =  0  plane.  In  fact,  it  appears  that  the  shed 
vortex  reaches  the  ground,  is  then  surrounded  by  the  crossflow,  and  thus  is  trapped  at  the  ground. 
Subsequently,  this  vortex  ring  is  lifted  off  the  ground  in  the  outflow  direction,  as  can  be  seen  in  the 
three-dimensional  perspective  (Figure  24).  Visualization  at  various  other  times  in  the  forcing  cycle 
have  shown  that  the  picture  in  Figure  25c  does  not  change  significantly. 
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Figure  24.  The  three-dimensional  perspective  of  the  vorticity  for  normal  impinging  jet  with 
crossflow  (case  B)  at  t  =  13.1.  Jet  forced  at  Si  =  0.47. 
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i.  a  — % orticity  in  (he  y— r  plane  at  x  =  0  at  t  =  10.8 


b.  a  orticity  in  the  y—z  plane  at  x  =  0  at  t  =  12.2 


c.  y -vorticity  in  the  x—z  plane  at  y  =  0  at  t  =  12.2. 
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d.  y— vorticity  in  the  x—z  plane  at  y /£>  =  1.5  at  t  =  12.2 


Figure  25.  The  vorticity  contours  in  the  flow  field  for  the  forced  normal  impinging  jet  with 
crossflow  (case  B).  Contour  interval  is  0.5  for  figures  25a-25c  and  0.1  for  figure  25d. 
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The  v-vorticity  in  the  y/D  =  1.5  plane  (Figure  25d)  shows  complex  vortical  motion  occurring  at 
this  location.  The  two  vortical  structures  seen  in  this  figure  should  not  be  confused  with  the  two 

streamwise  vortex  tubes  seen  in  Figure  24.  This  is  because  the  full  jet  in  Figure  24  is  shown  by  Hip¬ 
ping  the  computed  C  <  ^  /D  <  1.5  domain  to  the  -1.5  <  x  /D  <  0  domain  by  using  the  symmetry  con¬ 
dition  at  x  —  0  In  the  present  figure  (Figure  25d),  only  one-half  of  the  full  jet  is  shown  (i.e.. 

0  <  x  /D  <  1.5  ).  Thus,  the  vortical  structure  on  the  left  in  this  figure  corresponds  to  the  streamw  ise 
ortex  tube  in  Figure  24,  and  the  structure  on  the  right  is  the  cross-sectional  image  of  the  curved  vor¬ 
tex  ring  from  the  front  of  the  jet.  The  kidney-shaped  structure  in  Figure  24b  would  be  formed  by  the 
vortex  to  the  right  of  the  x  =  0  axis  and  its  image  to  the  left  of  the  x  =  0  axis.  Furthermore,  due  to 
symmetry,  it  is  clear  that  the  image  will  have  positive  uy  vorticity,  which  is  opposite  to  the  negative  u-y 
vorticity  in  this  figure.  By  reviewing  this  figure  along  with  the  three-dimensional  perspectives  shown 
in  Figure  24,  these  features  can  be  visualized. 

Figure  26  shows  the  velocity  vector  fields  in  the  various  planes  for  this  simulation.  Figures  26a 
and  26b  are  at  the  same  times  as  shown  in  Figures  25a  and  25b.  The  presence  of  a  rotating  region 
ahead  of  the  jet  is  clearly  seen.  However,  through  comparison  with  the  vorticity  contours  at  the  same 
location  and  time,  it  is  clear  that  whereas  the  recirculation  region  is  quite  large,  the  associated  |w*i  vor¬ 
ticity  in  this  region  is  quite  small.  In  experimental  studies,  the  presence  of  a  vortex  is  usually  identi¬ 
fied  from  the  streakline  patterns.  This  corresponds  to  the  velocity  vector  field  as  shown  in  these  fig¬ 
ures  However,  as  mentioned  earlier,  to  identify  this  region  as  a  vortical  region  would  not  be  correct, 
as  can  be  seen  by  comparing  Figure  25  with  this  figure. 

In  the  plane  of  the  fountain  at  x  /D  =  1.5  (Figure  26c),  the  flow  is  mostly  in  the  outflow  direction 
and  thus  there  is  no  significant  upwash  occurring.  Although  there  is  some  upward  motion  in  the  How 
field,  this  motion  is  turned  back  to  the  outflow  direction  as  the  outflow  boundary  is  approached. 

Figures  26d  and  26e  show  the  velocity  vector  field  in  the  x-z  plane  at  y  =  0  and  y  /D  =  1.5, 
respectively,  and  correspond  to  the  vorticity  fields  shown  in  Figures  25c  and  25d.  It  can  be  seen  that 
the  impinging  jet  starts  to  spread  toward  the  fountain  but,  due  to  the  effect  of  the  crossfiow,  the  velo¬ 
city  vectors  change  direction  and  are  predominantly  in  the  y-direction,  which  is  out  of  the  frame  of 
this  view.  The  recirculation  caused  by  the  streamwise  vorticity  behind  the  jet  and  the  vortex  ring 
from  the  ground  in  the  x-z  plane  can  be  seen  in  Figure  26e  quite  clearly.  Again,  the  recirculation  on 
the  left  near  the  x  =  0  axis  corresponds  to  the  streamwise  vortex  tube  related  to  the  kidney-shaped 
structures  observed  in  the  vorticity  plots;  the  motion  on  the  right  is  due  to  the  vortex  ring  lifting  off 
the  ground.  Furthermore,  due  to  symmetry  at  x  =  0,  the  recirculation  to  the  left  of  the  x  =  0  line 
would  be  opposite  to  the  recirculation  observed  to  the  right  of  the  x  =  0  axis;  thus,  the  recirculation 
pattern  associated  with  the  streamwise  vortex  tubes  is  counterrotating.  This  is  again  in  good  agree¬ 
ment  with  experimental  observations.4,22 

Figures  27a  and  27b  show  the  pressure  contours  in  the  y-z  plane  at  x  =  0  at  the  same  times  as  arc 
shown  in  Figures  25  and  26.  The  peak  stagnation  pressure  occurs  slightly  downstream  of  the  jet 
centerline  (i.e.,  y  =  0)  at  a  location  0.2 ID.  This  is  quite  similar  to  the  value  obtained  even  when  there 
is  a  fountain  present  (Figure  20c).  The  peak  value  and  the  location  are  not  significantly  affected  by 
the  forcing.  The  pressure  contours  in  the  x-y  plane  at  z  /H  =  0  (i.e.,  the  ground  plane)  are  shown  in 
Figure  27c.  The  compression  in  the  front  of  the  jet  due  to  crossfiow  can  be  seen.  This  is  again  as 
expected  from  previous  experimental  observations. 
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Figure  26.  The  velocity  vector  fields  in  the  flow  field  for  the  forced  normal  impinging  jet  with 
crossflow  (cont.). 
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b.  v  — r  plane  at  .v  =  0  at  t  =  12.2. 
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Figure  27.  The  pressure  contours  in  the  flow 
field  for  the  forced  normal  impinging  jet  with 
crossflow  (case  B).  Contour  internal  is  0.02. 


c.  ,v  —v  plane  at  r  /H  =  0.0  at  t  =  12.2. 
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Typical  velocity  profiles  in  the  flow  field  are  shown  in  Figure  28.  Figure  28a  shows  the  v-velocit\ 
profiles  in  the  y-z  plane  at  x  =  0  for  the  two  times  shown  earlier.  The  incoming  crossflow  velocity  is 
unaffected  by  the  effect  of  forcing  of  the  impinging  jet,  as  was  expected.  At  around  y  /D  =  0.33. 
there  is  definite  negative  v-velocity  near  the  wall.  This  is  related  to  the  recirculation  region  seen  in 
the  velocity  vector  diagrams.  In  the  positive  y-direction,  the  wall  jet  growth  is  similar  to  the  growth 
observed  in  the  earlier  simulation  (cf.  Section  5.3.1). 

Figure  28b  shows  the  corresponding  w-velocity  profiles  in  the  x-z  plane  at  y  =  0.  The  wall  jet 
grows  toward  the  fountain  axis,  but  due  to  the  crossflow  effect,  the  peak  values  are  quite  small.  There 
is  not  much  evidence  of  the  reverse  wall  jet  on  the  top  wall  near  x  /D  =  1.17,  which  is  close  to  the 
fountain  axis  at  x  /D  =  1.5. 

Figure  28c  shows  the  w-velocity  profiles  in  the  impinging  jet  for  the  two  times  of  the  forcing  cycle. 
The  effect  of  forcing  does  not  significantly  affect  the  jet  velocity  profile  except  near  the  potential 
core.  The  spreading  of  the  jet  is  quite  clear  and  is  similar  to  that  observed  earlier.  In  the  upstream 
direction  (y  /D  <  0),  there  are  hardly  any  noticeable  differences  in  the  velocity.  In  the  downstream 
direction  (v  /D  >  0),  there  is  some  positive  w-velocity  near  the  ground  which  may  be  related  to  the 
lift-off  of  the  vortex  ring,  as  seen  in  the  three-dimensional  perspective. 

Finally,  Figure  29  compares  the  wall  pressure  distribution  in  the  y-z  plane  at  x  =  0  for  the  forcing 
case  to  that  for  the  normal  impingement  case  (Section  5.2).  The  peak  pressure  at  the  stagnation  point 
for  the  case  with  crossflow  is  shifted  downstream  of  the  normal  impingement  case.  The  peak  values 
arc  slightly  affected  by  the  forcing.  Near  y  /D  =  -0.2,  there  is  an  increase  in  the  pressure.  Com¬ 
parison  with  the  velocity  vector  figures  indicates  that  this  increase  is  related  to  the  recirculation 
region  seen  ahead  of  the  impinging  jet.  The  pressure  is  higher  in  the  upstream  region  {y  /D  <  0)  when 
compared  to  the  downstream  region  (y/D  >  0).  This  is  partly  due  to  the  slowdown  of  the  crossflow  as 
the  jet  is  reached,  resulting  in  an  adverse  pressure  gradient  in  the  region  y  /D  <  0,  whereas  down¬ 
stream  of  the  jet  impingement  region,  the  crossflow  velocity  component  is  added  to  the  wall  jet. 

5.4  The  Normal  Impinging  Jet  in  Crossflow  with  Moving  Ground  Plane 

A  simulation  was  also  carried  out  for  the  case  when  the  ground  plane  is  also  moving  at  the  same 
speed  as  the  incoming  crossflow.  As  mentioned  before,  this  is  more  representative  of  a  V/STOL  air¬ 
craft  in  forward  motion.  The  cases  discussed  in  Section  5.3  correspond  more  closely  to  the  case  of  an 
experimental  study  in  a  wind  tunnel  where  the  ground  plane  is  fixed.  As  reported  in  the  literature, 
there  are  some  significant  differences  between  the  two  cases  that  must  be  taken  into  account  before 
the  wind  tunnel  data  can  be  utilized  to  describe  the  effects  of  crossflow  in  actual  flight  situations. 

The  simulation  performed  in  this  case  paralleled  the  case  discussed  in  Section  5.3.1,  in  which  the 
crossflow  was  begun  after  the  jet  had  impinged  on  the  ground.  The  results  were  acquired  at  the  same 
time  as  in  Section  5.3.1  so  that  a  direct  comparison  of  the  differences  in  the  flow  field  can  be  made. 

The  major  difference  is  expected  to  be  near  the  wall  region  since  in  the  present  case  the  ground 
moves  with  the  same  velocity  as  the  crossflow  and  thus  no  wall  boundary  layer  is  present.  We  must 
point  out  that  only  preliminary  conclusions  can  be  derived  from  this  study  since  the  low  grid  resolu¬ 
tion  near  the  wall  precludes  detailed  analysis  of  the  flow  there.  High-resolution  simulations  must  be 
performed  for  the  cases  described  in  Section  5.3.1  and  in  this  section  before  the  details  of  the  differ¬ 
ences  can  be  understood.  As  shown  below,  however,  even  with  low-resolution  computations,  signifi¬ 
cant  differences  arc  observed  that  are  in  remarkable  agreement  with  experimental  observations. 
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c.  »v -velocity  profiles  in  the  y-z  plane  at  x  =  0. 

Figure  28.  The  velocity  profiles  in  the  flow  field  for  the  forced  normal  impinging  jet  with 
crossflow  (cont.).  - t  =  10.8  -  •  i  =  12.2 


Figure  29.  The  ground  pressure  variation  for  the  forced  normal  impinging  jet  w  ith  crossflo 
(case  B).  □  /  =  9.4  Of  =  10.8  A  f  =  12.2 
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Figure  30a  shows  the  three-dimensional  vorticity  field  for  the  case  with  a  moving  wall  at  the  same 
instant  shown  in  Figure  18a  and  using  the  same  vorticity  level  and  viewing  angle.  Figures  30b  and 
30c  show  the  IwJ  and  components,  respectively,  for  the  moving  wall  case.  Comparison  with  the 
earlier  simulation  without  the  moving  wall  (Figure  18)  reveals  significant  differences.  Note  that  the 
"star-shaped"  vorticity  on  the  top  is  related  to  the  vorticity  on  the  top  wall.  As  expected,  there  is  very 
little  vorticity  present  in  the  near  wall  region  due  to  the  wall  motion.  The  vortex  ring  in  the  wall  jet 
seen  in  the  earlier  simulation  without  a  moving  wall  is  no  longer  present.  The  spreading  of  the  imp¬ 
inging  jet  shear  layer  now  shows  a  two-lobed  structure.  The  character  of  the  fountain  is  also  changed 
significantly.  For  example,  the  component  shows  that  the  structures  in  the  fountain  are  pushed 
farther  in  the  positive  y-direction  and  the  vortical  regions  in  the  fountain  do  not  extend  to  the  ground. 
This  is  understandable,  since  there  is  no  viscous  boundary  layer  on  the  ground  and  thus  the  "inviscid" 
wall  layer  acts  as  a  boundary  for  the  vorticity  in  the  flow  field. 

Figure  31a  shows  the  x-vorticity  in  the  y-z  plane  at  x  =  0.  Figure  31b  shows  the  corresponding 
velocity  vector  field.  Comparison  with  the  earlier  simulation  without  the  moving  wall  (Figures  19a 
and  20a,  respectively)  shows  that  the  effect  of  the  moving  wall  is  to  reduce  the  vortical  structure 
ahead  of  the  jet.  The  jet  shear  layer  bends  by  approximately  the  same  amount  even  with  the  moving 
wall.  This  is  understandable  because  the  dynamics  of  the  impinging  jet  should  not  be  significantly 
affected  by  the  moving  wall  except  in  the  near  wall  region.  Furthermore,  there  is  no  vorticity  in  the 
wall  region  due  to  the  absence  of  the  wall  boundary  layer.  The  velocity  vector  field  clearly  shows  that 
there  is  no  wall  boundary  layer  present  and  that  the  recirculating  region  ahead  of  the  impinging  jet  is 
reduced  in  size  considerably.  This  is  exactly  the  observation  made  in  the  experimental  studies,  where 
it  was  shown  that,  in  the  symmetry  plane,  the  vortex  ahead  of  the  impinging  jet  is  reduced  by  about 
30  percent  due  to  the  moving  wall.28  The  present  computations  indicate  a  good  agreement  with  this 
experimental  observation.  Figure  31c  shows  the  velocity  vector  field  in  the  y-z  plane  at  x  /D  =  1.5. 
There  is  clearly  an  inviscid  layer  present  near  the  wall,  in  contrast  to  the  wall  boundary  layer  seen 
when  the  wall  was  not  moving.  This  inviscid  layer  acts  as  a  boundary  in  the  plane  of  the  fountain 
and,  as  mentioned  earlier,  does  not  allow  the  fountain  to  evolve  in  the  near  wall  region. 

Figures  32a  and  32b  show  the  y-vorticity  contours  in  the  x-z  plane  at  y  =  0  and  y  /D  =  1.5, 
respectively.  In  comparing  these  figures  with  the  equivalent  figures  without  the  moving  wall  (Figures 
19b  and  19c),  it  is  clear  that  there  is  no  vortical  structure  in  the  fountain  axis  at  y  =  0.  This  is  due  to 
the  fact  that  the  fountain  is  pushed  farther  downstream  due  to  the  motion  of  the  inviscid  layer.  This 
was  also  seen  in  Figure  30.  The  vortical  structures  in  the  y  /D  =  1.5  plane  are  more  complex  with  the 
moving  wall.  The  core  of  the  vortex  in  the  fountain  is  now  closer  to  the  wall,  as  seen  on  the  right  side 
of  this  figure.  The  structure  on  the  left  is  the  image  of  the  vortical  layer  behind  the  impinging  jet  and 
is  associated  with  the  spreading  observed  in  Figure  30.  Note  that  there  is  some  positive  uy  vorticity 
between  the  two  structures.  This  is  possibly  due  to  the  wall  inviscid  layer  being  engulfed  by  the  vorti¬ 
cal  motions  in  the  fountain  region  and  the  jet  shear  region. 

Figures  32c  and  32d  show  the  corresponding  velocity  vector  fields  in  the  y-z  plane  at  y  =  0  and 
y  /D  =  1.5,  respectively.  The  velocity  field  in  the  y  =  0  plane  shows  negligible  upward  motion  in  the 
fountain,  and  there  are  no  significant  indications  of  the  wall  jet  near  the  ground  wall.  This  is  again  a 
result  of  the  wall  motion  in  the  y-direction,  which  causes  the  flow  field  in  the  near  wall  region  to 
move  primarily  in  the  positive  y-direction.  In  the  y  /D  =  1.5  plane,  the  recirculation  due  to  the  vorti¬ 
cal  motion  seen  in  Figure  32b  can  be  seen.  There  is  some  fluid  motion  to  the  fountain  region  (i.e., 
x  /D  =  1.5)  near  the  wall,  but  it  is  relatively  quite  weak. 
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a.  3D|w|  at  a  level  of  1.3 


b.  301^*  at  a  level  of  1.0 


c.  Velocity  vector  field  at  x  /D  =  1.5. 


Figure  31.  The  vorticity  and  velocity  fields  in  the  y—z  plane  for  the  normal  impinging  jet  with 
crossflow  and  w  ith  the  moving  ground  plane  at  t  =  13.1. 
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b.  3»— vorticity  at  y  /D  =  1.5.  Contour  interval  is  0.2. 


Figure  32.  The  vorticity  and  velocity  fields  in  the  x—  z  plane  for  the  normal  impinging  jet  with 
crossflow  and  with  the  moving  ground  plane  at  t  =  13.1. 
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Figure  32.  The  vorticity  and  velocity  fields  in  the  x-z  plane  for  the  normal  impinging  jet  with 
crossflow  and  with  the  moving  ground  plane  at  t  *  13.1  (cont.). 


TR-470/02-89 


55 


Figures  33a  and  33b  show  the  pressure  contours  on  the  top  and  bottom  walls,  respectively.  There 
is  a  low-pressure  region  on  the  top  wall  in  the  jet  shear  layer,  and  the  effect  of  crossflow  is  to  spread 
the  pressure  contours.  Note  that  the  pressure  increases  in  the  plane  of  the  fountain  but  the  high  pres¬ 
sure  is  shifted  further  toward  the  outflow  ( y/D  >  2).  This  is  the  region  where  the  fountain  reaches 
the  top  wall,  and  the  shift  downstream  is  due  to  the  crossflow.  The  pressure  contours  on  the  bottom 
wall  show  that  the  stagnation  pressure  occurs  downstream  of  the  y  =  0  location  by  about  0.28D,  and 
again  the  compression  effect  of  the  crossflow  can  be  seen.  There  is  a  low-pressure  region  behind  the 
jet. 

Figure  34  shows  typical  velocity  profiles  in  the  flow  field.  The  velocity  profiles  for  the  case  with  a 
nonmoving  wall  (Section  5.3.1)  at  the  same  time  is  also  shown  in  these  figures  to  allow  a  direct  com¬ 
parison  and  interpretation  of  the  effect  of  a  moving  wall.  Figure  34a  shows  the  v-velocity  profiles  in 
the  r-i  plane  at  .v  =  0.  The  lack  of  a  wall  boundary  layer  for  the  moving  wall  case  can  be  clearly  seen 
in  this  figure.  Further  away  from  the  wall,  the  motion  in  the  positive  y-direction  is  not  affected  signi¬ 
ficantly  by  the  moving  wall.  Figure  34b  shows  the  z<-velocity  profiles  in  the  x-z  plane  at  y  =  0.  As 
mentioned  earlier,  the  moving  wall  causes  significant  motion  toward  the  outflow  in  the  y-direction. 
Thus,  the  wall  jet  region  in  the  direction  of  the  fountain  axis  is  quite  weak.  Another  result  is  that 
there  is  no  significant  upwash,  which  in  turn  results  in  a  reduced  reverse  flow  on  the  top  wall. 

Figure  34c  shows  the  w-velocity  profiles  in  the  y-z  plane  at  x  —  0.  The  impinging  jet  is  not 
affected  significantly  by  the  moving  wall,  and  the  sweepback  of  the  jet  due  to  crossflow  occurs  in 
nearly  the  same  manner  as  without  a  moving  wall.  Ahead  of  the  jet,  however,  there  was  significant 
upwash  (positive  w-velocity)  for  the  nonmoving  wall  case,  which  was  shown  earlier  to  be  associated 
with  the  vortex  ancad  of  die  jet.  The  moving  wall  case,  on  the  other  hand,  shows  much  smaller 
upward  motion,  and  this  upwash  is  also  located  closer  to  the  jet  itself.  This  reduced  motion  is  related 
to  the  smaller  vortical  structure  and  recirculation  seen  in  earlier  figures  and  is  consistent  with  experi¬ 
mental  observations. 

Finally,  Figure  35  shows  the  wall  pressure  distribution  in  the  y-z  plane  at  x  -  0  for  three  cases. 
The  normal  impingement  case  and  the  nonmoving  wall  case  are  compared  in  this  figure  with  the  mov¬ 
ing  wall  case.  The  location  of  the  peak  stagnation  pressure  on  the  wall  with  crossflow  occurs  at  almost 
the  same  location  with  and  without  the  moving  wall.  In  general,  the  pressure  variation  on  the  wall  is 
not  significantly  affected  by  the  moving  wall. 
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a.  Top  wall. 


b.  Ground  plane. 


Figure  33.  The  pressure  contours  in  the  x-y  plane  for  the  normal  impinging  jet  with  crossflow 
and  with  the  moving  ground  plane  at  t  «  13.1.  Contour  interval  is  0.02. 
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6.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  report  presents  the  results  of  large-eddy  simulations  of  various  flow  fields  associated  with 
impinging  jets.  The  results  clearly  indicate  that  the  primary  technical  objectives  of  the  Phase  I  stud> 
have  been  successfully  achieved.  The  original  LES  code  for  studying  normal  impinging  jets  was  suc¬ 
cessfully  modified  to  simulate  inclined  jets,  normal  impinging  jets  in  crossflow,  and  normal  impinging 
jets  in  crossflow  with  a  moving  ground  plane.  Thus,  we  now  have  a  simulation  capability  that  can  be 
used  to  study  the  detailed  dynamics  of  fluid  motion  associated  with  various  flight  configurations  of  a 
V  STOL  aircraft  in  ground  effect.  These  modifications  were  accomplished  in  such  a  manner  that  the 
combination  of  various  flight  situations  can  also  be  simulated.  For  example,  a  V/STOL  aircraft  in 
pitch  while  in  transition  from  hover  to  forward  motion  can  also  be  simulated. 

Other  interesting  situations  that  can  be  simulated  are  the  effects  of  forward  acceleration  while 
moving  from  hover  tn  climb,  and  of  forward  deceleration  while  in  transition  from  forward  motion  to 
hover.  Another  possible  configuration  is  an  aircraft  undergoing  an  unsteady  pitching  motion.  The 
modifications  carried  out  in  Phase  I  were  made  to  take  these  features  into  account.  Although  these 
fight  situations  were  not  simulated  in  this  Phase  I  study,  the  success  of  the  present  computations 
leaves  no  doubt  of  the  possibility  of  studying  now  fields  associated  with  such  complex  night  confi¬ 
gurations.  It  is  worth  pointing  out  here,  however,  that  these  highly  unsteady  V/STOL  motions  occur 
at  very  high  frequency,  whereas  the  large-scale  vortical  motion  occurs  at  much  lower  frcqimn<-\ 
Therefore,  parametric  variation  of  the  impingement  angle  or  the  forward  speed  will  probably  be  suffi¬ 
cient  to  understand  the  complex  dynamics  of  the  fiuid  fiow.  Thus,  the  LES  capability  demonstrated 
in  this  Phase  I  study  may  provide  a  new  numerical  "experimental"  tool  that  can  be  utilized  to  under¬ 
stand  the  dynamics  of  the  flow  field  in  a  systematic  manner. 

The  basic  numerical  scheme  was  used  in  earlier  studies  to  investigate  the  effect  of  coherent  and 
random  excitation  of  normal  impinging  jets.  The  results  were  shown  to  agree  quite  well  with  experi¬ 
mental  observations  and  thus  the  validity  of  the  basic  code  was  demonstrated.  In  the  present  Phase  1 
study,  only  coarse  mesh  calculations  were  possible  due  to  resource  limitations.  However,  the  grid  was 
chosen  so  that  some  comparison  with  experimental  data  could  be  accomplished.  The  results  presented 
in  this  report  clearly  demonstrate  that  this  LES  code  is  capable  of  reproducing  features  in  the  How 
field  that  are  in  good  qualitative  agreement  with  some  of  the  experimental  observations. 

The  approach  taken  in  this  study  was  to  investigate  the  motion  of  large-scale  structures  in  the  flow 
field  by  exciting  the  impinging  jet  at  a  prescribed  frequency  in  the  jet  preferred  range.  This  ensures 
that  the  large  structures  generated  in  the  impinging  jet  will  be  controlled  by  the  forcing  frequency  and 
that  the  natural  instability  of  the  flow  will  be  suppressed,  resulting  in  a  "clean"  numerical  experiment. 

The  simulation  of  the  inclined  impinging  jet  showed  some  good  agreement  with  experimental 
observations.  The  oblique  impingement  causes  a  significant  portion  of  the  jet  fiuid  to  be  transported 
in  the  direction  of  inclination.  Theoretical  estimates  indicate  that  about  80  percent  of  the  mass  is 
transported  in  the  forward  direction.15,17  The  jet  impingement  point  is  behind  the  jet  center  line  for 
the  oblique  impingement  case,  which  is  in  agreement  with  experimental  observations.  In  addition,  the 
shift  forward  from  the  center  (y  =  0)  location  for  the  normal  impingement  case  was  also  determined  to 
agree  closely  with  the  experimental  value.  It  was  also  demonstrated  that  the  stagnation  streamline  for 
the  oblique  impingement  case  has  nonzero  vorticity  in  contrast  to  the  nonrotational  stagnation  stream¬ 
line  for  the  normal  impingement  case.  This  agreed  with  past  experimental  and  theoretical  observa¬ 
tions.17,18 
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The  effect  of  forcing  is  that  large  coherent  vortical  rings  arc  shed  at  the  forcing  frequency.  The 
phase  speed  of  these  structures  in  the  jet  was  estimated  to  be  around  0.6 U,  which  again  agrees  with 
experimental  values.3  The  presence  of  the  fountain  formed  by  the  collision  of  the  wall  jets  changes 
the  dynamics  of  the  structures  in  the  wall  region,  as  was  observed  earlier  in  the  normal  impinging  jet 
studies.  The  fountain  is  skewed  in  the  direction  of  inclination,  which  would  result  in  a  nonsymmctric 
upwash  and  consequently  an  unequal  lift-off  effect  on  the  top  wall. 

The  effect  of  crossflow  with  and  without  a  moving  wall  was  also  studied,  and  the  results  are  again 
in  agreement  with  some  experimental  observations.  For  example,  the  experimentally  observed4  for¬ 
mation  of  counterrotating  structures  behind  the  impinging  jet  was  also  numerically  computed.  How¬ 
ever,  it  was  shown  that  the  effect  of  the  proximity  of  the  fountain  (which  is  absent  in  most  of  the 
experimental  studies)  can  modify  the  dynamics  of  the  vortical  motion  significantly.  The  streamline 
compression  ahead  of  the  jet  due  to  crossflow  was  also  demonstrated  and  was  in  qualitative  agreement 
with  experimental  observations. 

The  impinging  jet  acts  as  an  obstacle  to  the  crossflow,  causing  the  crossflow  to  go  around  it.  The 
ring  vortex  ahead  of  the  impinging  jet  is  pushed  back  and  significant  three-dimensional  vortex  stretch¬ 
ing  occurs  due  to  crossflow.  The  total  vorticity  in  the  flow  field  is  redistributed  ir.  r  very  complex 
manner,  as  was  also  determined  in  experimental  studies.4,41  However,  this  redistribution  is  quite  dif¬ 
ferent  if  the  fountain  is  formed  (Section  5.3.1)  as  compared  to  the  case  when  the  crossflow  is  initiated 
before  tne  fountain  is  formed  (Section  5.3.2).  More  detailed  and  long-time  simulations  are  necessary 
to  evaluate  the  dynamics  of  the  large-scale  motion  under  crossflow.  These  aspects  will  be  considered 
in  more  detail  in  the  next  phase. 

The  effect  of  allowing  the  ground  plane  to  move  is  that  the  wall  boundary  layer  is  absent  and  an 
in  viscid  layer  is  present  close  to  the  wall.  This  results  in  the  reduction  of  the  vortex  size  and  strength 
ahead  of  the  impinging  jet,  which  again  agrees  with  experimental  observations.27,28  The  computations 
of  impinging  jets  in  crossflow  with  a  moving  ground  plane  presented  in  this  study  appear  to  bo  the 
first  attempt  to  numerically  simulate  these  flow  fields.  The  preliminary  results  obtained  indicate  good 
qualitative  agreement  with  experimental  observations. 

In  conclusion,  the  LES  capability  developed  and  demonstrated  in  this  Phase  I  study  clearly  shows 
that  all  the  technical  objectives  of  the  study  have  been  successfully  achieved.  The  numerical  code  can 
now  simulate  a  variety  of  V/STOL  flight  situations  in  ground  effect.  These  include  pitch  up  or  down, 
forward  motion  at  a  steady  velocity  and  under  forward  acceleration/deceleration.  The  effect  of  the 
moving  ground  plane  can  also  be  modeled  by  the  LES  code.  The  jet-to-ground  distance  and  the 
lateral  spacing  between  the  impinging  jets  can  also  be  varied  to  study  the  effects  of  lift-off  and  the 
dynamics  of  the  fountain.  Most  of  these  configurations  can  also  be  combined  to  study  realistic 
V/STOL  flight  in  ground  effect.  For  example,  an  actual  flight  profile  of  the  V/STOL  near  t he  ground 
may  be  composed  of  a  lift-ofT  from  the  ground,  a  pitch-up  followed  by  a  rapid  acceleration  to  a  con¬ 
stant  velocity  with  which  the  aircraft  would  climb  out  of  ground  effect.  The  LES  code  demonstrated 
in  this  Phase  I  study  can  be  used  to  study  most  of  this  flight  profile  in  a  spatially  and  temporally  accu¬ 
rate  manner.  The  dynamics  of  large-scale  vortical  structures  in  the  impinging  jets,  the  wall  jets  and  in 
the  fountain  can  be  evaluated  as  demonstrated  here. 

Further  simulations  using  high  grid  resolution  are  required  to  understand  the  various  features  of 
the  complex  flow  field  associated  with  V/STOL  flows,  especially  at  high  Reynolds  numbers.  The 
effect  of  varying  the  jet-to-ground  distance,  H,  and  the  spacing  between  the  impinging  jets  must  also 
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be  evaluated  since  these  parameters  control  the  dynamics  of  the  fountain,  which  in  turn  controls  the 
lift-off  effect  on  the  V/STOL  aircraft.  Another  important  aspect  of  V/STOL  aircraft  in  ground  effect 
is  the  entrainment  of  dust  from  the  ground,  which  can  be  reingested  into  the  engine  intakes  and  possi¬ 
bly  cause  engine  malfunctions.  In  principle,  by  linking  a  "particle  tracking"  algorithm  to  the  LES 
code,  it  would  be  possible  to  study  the  motion  of  "particles"  injected  from  the  ground  in  order  to 
better  understand  the  entrainment  mechanism.  These  aspects  can  be  considered  in  the  next  phase. 

Finally,  it  must  be  emphasized  that  the  large-eddy  simulation  technique  provides  all  the  flow  vari¬ 
ables  of  interest,  including  properties  such  as  vorticity,  which  is  difficult,  if  not  impossible,  to  measure 
experimentally.  In  principle,  with  proper  data  management,  these  properties  can  be  recorded  in  the 
entire  computational  domain  for  the  entire  time  length  of  the  simulation.  Thus,  additional  information 
such  as  the  turbulent  kinetic  energy  and  Reynolds  stress  variations  in  the  flow  field  can  be  computed, 
as  demonstrated  in  an  earlier  study.10  The  temporal  record  of  the  velocity  and  pressure  fluctuations 
in  the  flow  can  also  be  used  for  spectral  analysis,  which  can  provide  additional  information  on  the  fre¬ 
quency  content  of  the  fluctuations.  Detailed  analysis  of  the  computed  results  and  comparison  with 
experimental  data  can  then  be  carried  out,  as  will  be  proposed  for  the  next  phase. 
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